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I. Introducere 


I.1. Etape în utilizarea programelor profesionale de analiza cu elemente finite 


Etapele principale în utilizarea programelor profesionale de analiză cu elemente finite a 
unei probleme sunt reprezentate schemetizat în Fig. 1. 


Definirea Analiza, 
problemei Decizie 


Preprocesare Procesare Postprocesare 


Fig. 1 


Preprocesare: 
e Generarea geometriei modelului; 


e Stabilirea proprietăților de material şi a caracteristicilor geometrice; 
e Discretizarea modelului cu tipul de elemente finite adecvate problemei; 
e Aplicarea condițiilor la limită pentru modelul cu elemente finite; 
e Definirea încărcărilor pe model. 
Procesare: 


e Verificarea datelor de intrare; 
e Lansarea în execuţie a modelului complet. 


Postprocesare: 
e Reprezentarea rezultatelor (deplasări, tensiuni, deformatii, etc.) sub formă de grafice, 
liste, contururi (isocurbe); 
e Analiza şi interpretarea rezultatelor. 


I.2. Prezentarea pachetului software de analizá cu elemente finite COSMOS/M 


COSMOS/M este un sistem complet, modular de calcul cu elemente finite, dezvoltat de 
firma Structural Research and Analysis Corporation pentru statii de lucru si calculatoare PC. 
Este destinat modelárii unei game largi de probleme prin metoda elementelor finite, capabil sá 
rezolve probleme liniare şi neliniare (statice şi dinamice), probleme de calcul la oboseală, 
probleme de optimizare, probleme de mecanica fluidelor, transfer de căldură şi electromagnetice. 

Modulul principal al programului, denumit Geostar, este destinat preprocesării cât şi 
postprocesării (generarea modelului geometric, discretizarea, aplicarea încărcărilor şi a 
condiţiilor la limită, prezentarea rezultatelor) prin intermediul lui realizându-se şi comunicarea 
cu utilizatorul. Programul mai conţine pe lângă acest modul şi o serie de module care controlează 
diferite categorii de analize (vezi figura). 

Geometria modelului poate fi creată cu ajutorul modulului Geostar, sau poate fi 
importată din diverse medii de dezvoltare a aplicaţiilor CAD sau din alte medii de calcul cu 
elemente finite. 


SolidWorks 

Microstation 

Modeler 

Solid/Edge 
CAD IGES 
CAD DXF 
NASTRAN 

Import date Export date ANSYS 


PATRAN 
ADAMS 
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Pro/ENGINEER 
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ANSYS 
PATRAN HSTAR 
(transfer de cáldurá) 


NSTAR 


(analizá neliniará staticá si dinamicá) 


STAR 
(analizá liniará staticá) FSTAR 
(calcul la oboseală) 


DSTAR 
(analiză dinamică — frecvenţă, flambaj) a ASTAR 3 
(analiză dinamică avansată) 
Module FFE 
(analiză statică, termică, de frecvență) Der OP TSTAR Nr 
optimizare si senzitivitate 
(opt § tivitate) 


COSMOS/HFS 2D 
COSMOS/HFS 3D 


Module de baza 
Module optionale 


COSMOS/CAVITY 


FLOWPLUS 
(analiza curgerii fluidelor) 


Ecranul de lucru al modulului Geostar este format din trei ferestre distincte, şi anume: 
1. Fereastra principală (Main) în care se construieşte modelul (fereastra grafică) şi în a cărei 
parte superioară se află bara cu meniurile "pull-down": 


[C Geostar 2.95 (256k d 
Im Edit Geometr IsBC Control Display Analysis Results Windows Help 


FILE —grupează comenzi pentru operaţii cu fişierele programului şi imaginile captate. 

EDIT — comenzi pentru operaţii de identificare, renumerotare, listare, ştergere, etc. 

GEOMETRY - meniul contine comenzi de construire a modelului geometric. Se pot genera 
entitáti geometrice de tipul: puncte, curbe, suprafete, volume, contururi, regiuni, poliedre, párti. 
Pentru construirea geometriei modelului utilizatorul poate folosi mai multe sisteme de 
coordonate (carteziene, cilindrice, sferice). 

MESHING — grupează comenzile de discretizare în elemente finite a modelului. Cu ajutorul 
acestor comenzi se defineşte densitatea reţelei de elemente finite, se realizează discretizarea 
parametrică sau automată a diferitelor entităţi geometrice, se operează asupra nodurilor sau 
elementelor finite. 

PROPSETS — comenzile grupate în meniul PropSets asigură o parte importantă a preprocesării. 
Cu ajutorul lor se aleg tipul de elemente finite cu care se va lucra, se definesc proprietățile de 
material, se stabilesc constantele reale asociate elementelor (arie, momente de inerție, grosime, 
etc.). 

LOADSBC - grupează comenzile ce permit definirea încărcărilor şi condiţiilor la limită. 
Deplasările şi deplasările sunt precizate direct nodurilor individuale sau indirect prin intermediul 
entităților geometrice de care nodurile sunt legate. 

CONTROL — grupează comenzi ce asigură comanda şi controlul a numeroase activități în timpul 
şi după crearea modelului (controlul mediului grafic, activarea şi selectarea prin diferite 
modalităţi a anumitor seturi de entități, comunicarea cu sistemele CAD, operaţiile de output). 
DISPLAY — meniul conține comenzile pentru controlul afişării. 

ANALYSIS — cu ajutorul comenzilor din acest meniu se defineşte tipul analizei pentru problema 
solutionatá şi se setează opţiunile legate de rezolvarea ei. 

RESULTS — grupează comenzile ce asigură postprocesarea, vizualizarea-preluarea şi analiza 
grafică a rezultatelor. 

WINDOWS — comenzile pentru gestionarea ferestrelor. 

HELP — help-ul interactiv de lucru. 


=|) x] 2. Fereastra GeoPanel poate contine sau nu o parte din meniurile din 


Semel | File | fereastra Main şi în plus butoane pentru pentru comenzi de tip view, 
zoom, rotate, scale, status, care uşurează utilizatorului observarea şi 


. Meshing | ___Edt | codstiuitea modelului, 
PropSets | Control | 
LoadsBC | Display | 
Analysis | Results | 
Active Win [Main — v] 


New win | Tile | Cascade 


I ———A—— "EI 
Save/Restore... | 


DView | Demo Config 


Zoom Te 


EM HW 
cloa 
mi 
C x 


— E ž 
[v Console | Set Post | 
Status | Status2 | Status3 | 


Display menu 


3. Ín Geostar Console (fereastra text) este disponibilá o linie de comandá unde se pot introduce 
comenzi, direct sau prin citire din fişiere. Aici se acumulează progresiv lista de comenzi 
corespunzătoare acţiunilor de construire a modelului. În această fereastră programul mai afişează 
o serie de mesaje ajutătoare pentru utilizator, avertizându-l de exemplu dacă datele introduse 
sunt incorecte sau incomplete, anuntandu-l care este rezultatul anumitor acţiuni (numărul 
entităților legate sau renumerotate), etc. 


„oii 
NEW pu 
GeoStar > m 


a 


Programul poate fi comandat interactiv prin intermediul meniurilor, al butoanelor activate cu 
mouse-ul si al comenzilor introduse de la tastatură. Programul poate lucra şi prin încărcarea şi 
executarea unui fişier de comenzi de tipul *.geo sau *.ses, cu ajutorul comenzii Load file, fişier 
care poate fi editat off-line. 


Capitolul II 


ELEMENTE FINITE DE TIP BARÀ 


II.1. Elemente finite de tip bará cu un singur grad de libertate pe nod 
(elemente finite de tip TRUSS) 


II.1.1. Elemente finite de tip bará cu un singur grad de libertate pe nod in 
plan (elemente finite de tip TRUSS2D) 


Problema: Truss2D-1 

Caracteristicile de incárcare: forte aplicate in noduri 

Sá se determine deplasările nodale şi tensiunile care se produc în barele 
structurii, Fig. II.1.1, sub acţiunea forţei verticale F = 100 kN. Barele 1 - 6 au secțiune 
circulară cu diametrul d = 30 mm, iar barele 7 — 13 au secțiune dreptunghiulară cu 
dimensiunile bxh = 20x50 mm, şi sunt realizate din oţel cu modulul de elasticitate 
longitudinal E = 2:10? MPa. 


Fig. IL.1.1. Schema geometrică de încărcare şi rezemare a structurii plane din bare articulate 
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Rezolvarea problemei cu programul COSMOS/M 


Important! 
Datoritá volumului mare de date salvarea documentelor create prin deschiderea unei 
aplicatii noi trebuie sá se facá intr-un director de lucru separat, creat special pentru 
acest scop. 


Realizarea modelului (preprocesarea): 


1. Pregátirea pentru construirea modelului eres E 
Display» View Parameter? View | p și i E = i i i 
VIEW, 0, 0, 1, 0 Pim ici Ex ten ter Kom Keep Mark DMsh 
. LoT iv iv 
Se vor activa din: c vB E 0r rr 
Control? Utility» Status Table 1 sr [- m m = O m mj m 
E E VL CBI " m m m 

cásutele corespunzátoare pentru puncte, curbe, | rr EN 0r 
elemente şi noduri. rr Or a—dé 0 [- B 

PT Mr L3 of M 

PA C r BI 0 [- 

DEM BF E " - 

EL vri 0 [- 

o : 

Color Index 

E 2 3 4 5 6 7 8 9 1011 1213 14 15 16 

Ce CI  — — II] 

Help | Qut | 


2. Generarea punctelor (se definesc prin 
coordonatele lor) 


D aa Ata 
PT, 1, 0, 0,0 [pr x! 
PT, 2, 500, 0, 0 Keypoint [1 
PT, 3, 1000, 0, 0 X-Coordinate Value fo 
PT, 4, : ; " 2 i Y-Coordinate Value [0 
PT, 5, | 0, : [p] —— — — — 
PT. 6, 500, 800, 0 Z-Coordinate Value 
PT, 7, 1000, 800, 0 Help | _ Cancel | 


PT, 8, 1500, 800, 0 


Íncadrarea desenului in limitele ecranului se realizeazá cu ajutorul butonului AUTO 
din fereastra din stânga Geopanel, sau 

Display» Display_Option> Scale 

SCALE, 0 


x 

3. Generarea curbelor Curve E 

Geometry? Curves? Line with 2 Pts Keypoint1 [7 
CRLINE, 2, 2, 3 d 


CRLINE, 3, 3, 4 Hep | Cancel | 
CRLINE, 4, 4, 5 
CRLINE, 5, 6, 7 


E] 


4. Alegerea tipului de elemente finite 

Ín biblioteca de elemente a programului COSMOS elementul este denumit TRUSS2D 
(Thermal and Linear Structural 2D Spar/Truss Element). Este un element finit uniaxial 
cu două noduri, pentru analiza structurală şi termică a modelelor plane. 


E 
u group i 


Element Name [ANU o 2D truss/spar element 


[ Continue ] Help | Cancel | 
EM E 


DP1:Unused opionl[ — 
OP2:Unused option jf ©. 
OP3:Unused option jf ©. 
OP4:Unused option jf | 
OP5:Material Type [0: Linear Elastic >] 
OPE: Displacement Formulation [osma o 
OP?:Material Creep JO: Nocep v] 
OP8:Unused option jf ©. 


Help | Cancel | 


PropSets> Element Group 
EGROUP, 1, TRUSS2D, 0, 0, 0, 0, 0, 0, 0, 0 


5. Stabilirea proprietăților de material 
Pentru analiza structurală proprietatea de material ce trebuie indicată este modulul de 
elasticitate. 

PropSets> Material Property 

MPROP, 1, EX, 2e5 
Observaţie: Se poate apela şi la un material din biblioteca programului, dar setând cu 
atenţie unităţile de măsură. 
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Brickman x 
pu property set fi 


Material Name [A - STEEL: ALLOY STEEL -| 


Unit Label |FPS: British (Inch/Pound/Second) Y 


FPS: British (Inch/Pound/Second) 
MKS: Metric Centimeter/Kilogram/Second 


Sl... International (Meter/Kilogram/Second]) 


6. Stabilirea proprietátilor geometrice ale sectiunilor 

Din enuntul aplicației reiese cá structura analizată este formată din bare având două 
secţiuni diferite. Se vor defini în consecință două seturi de constante reale, câte unul 
pentru fiecare secțiune, acestea asociindu-se ulterior elementelor finite ce se vor 
genera. 

PropSets> Real Constant 

RCONST, 1, 1, 1, 1, 706.858 


[wr O NT 
Associated Element group il 


Real Constant set fi 
Help | Cancel | 


[]RcoNsT p 


Start location of the real constants J == 


[ Continue ] Help | Cancel | 


Liner es) 
Mo. of real constants to be entered E 


Help | Cancel | 


RCONST, 1, 2, 1, 1, 1000 L]RcONST x] 
RC1 : Cross-sectional area [706. 858 


Help | Cancel | 


7. Discretizarea barelor în elemente finite 

Discretizarea barelor structurii presupune două acțiuni, respectiv activarea seturilor de 
proprietăți (material şi geometria secțiunii) asociate barelor ce urmează a fi discretizate 
şi discretizarea propriu-zisă; în aplicația prezentată barele structurii sunt realizate din 
acelaşi material, dar au secţiuni diferite. 
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Activarea setului 1 de 


constante reale, asociat Set Label 8 FAC: Real Constant 


barelor 1 — 6: 
Control Activate» Set Entity Heb | NN 


ACTSET, RC, 1 


[J ACTSET xj 


Real constant set number il 
Help | Cancel | 


Discretizarea barelor 1 — 6: 


Meshing> Parametric Mesh» — mem xÍ 
Curves Beginning Curve [1 

M_CR, 4, 5, £2, 4,4 — | 
Activarea setului 2 de Increment [1 

constante reale, asociat barelor Number of nodes per element [2 z] 
7-13: Number of elements on each curve h 

Control? Activate? Set Entity Spacing ratio [1 

ACTSET, RC, 2 Keypoint to define principal axis if 3 nd/el | 


Discretizarea barelor 7 — 13: 
Meshing> Parametric_Mesh> 
Curves 

M CR, 7, 13, 1, 2, 1,1 


Hep | | Cancel | 


8. Legarea nodurilor suprapuse si 


renumerotarea lor 
Meshing> Nodes? Merge Beginning Node [1 
NMERGE, 1, 26, 1, 0.0001, 0, 1, 0 EndngNode [E —— ^ — —^ — 
Total number of nodes merged — 18 B 
Increment ooo 
Tolerance [o. 0001 
All/Among flag fo: Al >] 


[_]NCOMPRESS X| 
Edit» Compress> Nodes Beginning Node i 
NCOMPRESS, 1, 12 Ending Node [13 


Number of Nodes compressed — 4 
i P _ Hep | Cen | 
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9. Definirea condiţiilor la limită si a încărcării 


Ce ju i -—— E 


nodul re reprezint 
odul care da a Beginning Node Moo 


articulația plană din A 

LoadsBC>Stru Displacement label ene le X translation 7] 
ctural>Displacement> Vale 
Define by Nodes Ending Node [1 

DND, 1, UX, 0, Increment ————— ——— — — 
1, 1, UY; Additional Displacement labels if any (L2... LE] DY — 


Help | Cancel | 


Conditii la limitá in nodul care reprezintá reazemul mobil din din B 
LoadsBC> Structural? Displacement» Define by Nodes 


OS TT 0s 


Beginning Node 5 | 


Displacement label pea al 
Value [D ©. 
Ending Node fb © 
Increment n 
AddiionalDisplacementlabelsifany(L2...L6]| — 


Help | Cancel | 


usce ia "MEN x: 
LoadsBc»Structural^ Force» Define by PREIEI Îi 
Nodes 
FND, 2, FY, =] 00000, 2, 1; Force "E FY: Force in Y 

Value [1 00000 
Solutionarea problemei (procesarea): Ending Node [2 
10. Verificarea datelor de intrare Increment [1 


Analysis Run Check ln cec es 


R CHECK, STATIC Analysis type [Ef 
11. Lansarea in executie 


Analysis? Static» Run Static Analysis Help | Cancel | 


R STATIC 
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Prezentarea rezultatelor analizei (postprocesarea): 


12. Vizualizarea rezultatelor cu ajutorul comenzilor din lista Results 
Tensiunile normale o, din barele structurii: 


Resulis Plot» Stress "CI -——- x 
ACTSTR, 1, SX, 2, 1, 0, 0 Load case number We os 
Component [sx: Normal Stress (4 dir) >] 
Stress flag PR =] 


Layer number n | 5  — | 


Face flag (Shell) lc— Top =] 


Coordinate system [n 0 


[| Contour Plot ] Vector Plot | | sa Plat] Section Plat | Help | Cancel | 


Results? Plot» Stress? Contour Plot 

STRPLOT, 0, 1, 13, 1, 0, 176.407 E 

Line flag Jil 

Beginning Element m 

EndngElemnt[13 — 

Increment n 1 1 1 | 

Shape flag [0 Undetomed = 

Deformed scale factor [640 — — ™S 


Lin STRESS Lc-1 


Slgma_X 
100.000 
B 5 7 6 8 


76.4450 
52.8890 
29.3340 
5.77820 
-17.777 
-41.333 

-64.888 


-88.444 


Id 


Tensiunile normale Sigma x din barele structurii 
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Listarea deplasărilor din nodurile structurii: 
Results>List>Displacement/Response/Reaction 
DISLIST, 1, 1, 1, 8, 1, 0 


=) x) 


[]DISLIST,1,1,1,8,1,0 


1 
UX UY UZ RX RY RZ 
1 0.000e«000 0.000e+000 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
2 1.658e-001 -1.134e+000 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
3 3.316e-001 -6.043e-001 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
4 3.868e-001 -3.367e-001 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
5 4.421e-001 0.000e+000 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
6 3.868e-001 -7.337e-001 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
7 2.763e-001 -6.043e-001 0.000e+000 0.000e+000 0.000e+000 0.000e+000 
8 1.658e-001 -3.367e-001 0.000e+000 0.000e+000 0.000e+000 0.000e+000 


Forţele de reacțiune: 
Results>List>Displacement/Response/Reaction 
DISLIST, 1, 2, 1, 8, 1, 0 


[susti T 
Load case il 
Node RFX REY RFZ RFRES 


7 500e+004 0.000e+000 7.500e+004 
2.500e+004  0.000e+000 2.500e+004 
1.000e+005  0.000e«000 1.000e+005 


1 0.000e+000 
5 0.000e+000 
-O 000e+000 


NOTA: 

Datele referitoare la Preprocesare, Procesare si Postprocesare sunt salvate 
intr-o bazá de date constituitá dintr-un numár diferit de fisiere in functie de tipul 
analizei. 


In cazul problemei de fata iatá care sunt fisierele create: 


B Truss2D_1 


Truss2D 1.DIS [E] Truss2D_1.1LC E] Truss2D_1.wnt 


AË Truss2D_1 
E Truss2D 1 
|a) Truss2D_1 
|a) Truss2D_1 
El Truss2D 1.bc 
(9) Truss2D. 1.CH 


Truss2D_1.dsp 
Truss2D 1.egr 
Truss2D 1.ele 
Truss2D. 1.gcr 
Truss2D 1.gen 
Truss2D 1.gpt 


El Truss2D 1.csy Truss2D_1.IDA 


Truss2D. 1.LCD 
Truss2D. 1.LDS 
Truss2D 1.lod 

Truss2D 1.nod 
Truss2D 1.rlc 

Truss2D. 1.5TE 
Truss2D. 1.V50 


Fişiere de tip text: Truss2D 1.out — confine listare 


rezultate; Truss2D.ses — contine sesiune comenzi. 
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Fişierul de listare a sesiunii de comenzi (cu extensia: .ses) are următorul 
conținut în cazul aplicaţiei de fata: 


VIEW,0,0,1,0 
PLANE,Z,0,1 
PT,1,0,0,0 

PT,2,500,0,0 
PT,3,1000,0,0 
PT,4,1500,0,0 
PT,5,2000,0,0 
PT,6,500,800,0 
PT,7,1000,800,0 
PT,8,1500,800,0 
SCALE,0 
CRLINE,1,1,2 
CRLINE,2,2,3 
CRLINE,3,3,4 
CRLINE,4,4,5 
CRLINE,5,6,7 
CRLINE,6,7,8 
CRLINE,7,1,6 
CRLINE,8,2,6 
CRLINE,9,6,3 
CRLINE,10,3,7 
CRLINE, 1 1,3,8 
CRLINE, 12,4,8 
CRLINE, 13,8,5 
EGROUP, 1, TRUSS2D,0,0,0,0,0,0,0,0 
MPROP,1,EX,2E5 
RCONST, 1,1,1,1,706.856 
RCONST, 1,2,1,1,1000 
ACTSET,RC,1 

M CR,1,6,1,2,1,1 
ACTSET,RC,2 

M CR,7,13,1,2,1,1 
NMERGE,1,26,1,0.0001,0,0,0 
NCOMPRESS,1,12 
DND,1,UX,0,1,1,UY, 
DND,5,UY,0,5,1, 
FND,2,FY,-100000,2,1 
C* R CHECK,STATIC 
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Fişierul de listare a rezultatelor (cu extensia: .out) are următorul conținut in 
cazul aplicaţiei de fata: 


LILLLLLLLLLLLLLLLELLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLE 
LILLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLILLLLLLEL 
ek "ek 
2k "ek 


** CCCC OOOO SSSS M M OOOO SSSS /M M ** 
**C OOS MMMMOOS /MM MM ** 
** C OO SSSS MMMM OO SSSS / MMMM ** 
**C OO SMMMOO S/ MMM ** 

** CCCC OOOO SSSS M M OOOO SSSS/ M M ** 


dk ok 

ok ok 

NE VERSION: 2.5 Mt 

AE DISTRIBUTED BY: TX 

** STRUCTURAL RESEARCH AND ANALYSIS CORPORATION  ** 
ut 12121 WILSHIRE BLVD. SUITE 700 mum 

2n LOS ANGELES, CALIFORNIA 90025 cui 

AT TEL. NO. (310) 207-2800 kd 

id COPYRIGHT 1988 S.R. A.C. om 

dk ok 


LILLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLE 
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Problem name: Truss2D 1 

Date : 05/27/2009 Time: 14:01:29 
Title 

Subtitle 


CONTROL INFORMATION 
NUMBER OF LOAD CASES ........... (NLCASE) = 1 
SOLUTION MODE ............... (MODEX) - 0 
EQ. 0, STATIC ANALYSIS 
EQ. 1, BUCKLING ANALYSIS 
EQ. 2, DYNAMIC ANALYSIS 
SOLVER TYPE... nme (ISOL) = 0 
EQ. 0, DIRECT SPARSE SOLVER 
EQ. 1, DIRECT SKYLINE SOLVER 
EQ. 2, ITERATIVE SOLVER 


THERMAL LOADING FLAG............ (ITHERM) = 0 
EQ. 0, NO THERMAL EFFECTS CONSIDERED 
EQ. 1, ADD TEMPERATURE EFFECT 
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GRAVITY LOADING FLAG. ............ (IGRAV) = 0 
EQ. 0, NO GRAVITY LOADING CONSIDERED 
EQ. 1, ADD GRAVITY LOADING EFFECT 
CENTRIFUGAL LOADING FLAG. .......... (ICNTRF) = 0 
EQ. 0, NO CENTRIFUGAL LOADING CONSIDERED 
EQ. 1, ADD CENTRIFUGAL LOADING EFFECT 
IN-PLANE STIFFENING FLAG .......... (INPLN) = 0 
EQ. 0, NO IN-PLANE EFFECTS CONSIDERED 
EQ. 1, IN-PLANE EFFECTS CONSIDERED 
SOFT SPRING ADDITION FLAG......... (ISOFT) = 0 
EQ. 0, NO SOFT SPRING OPTION 
EQ. 1, SOFT SPRING ADDED 


SAVE DECOMPOSED STIFFNESS MATRIX FLAG ...(ISAVK) = 0 
EQ. 0, DO NOT SAVE DECOMPOSED K 
EQ. 1, SAVE DECOMPOSEDK 
FORM STIFFNESS MATRIX FLAG ......... (IFORMK) = 0 
EQ. 0, FORM STIFFNESS MATRIX 
EQ. 1, USE EXISTING DECOMPOSED STIFFNESS MATRIX 


SPIN SOFTENING FLAG............ (ISPIN) = 0 
EQ. 0, NO SPIN SOFTENING EFFECTS CONSIDERED 
EQ. 1, SPIN SOFTENING EFFECTS CONSIDERED 
INERTIA RELIEF FLAG ............ (IFORMK) = 0 
EQ. 0, NO INERITA RELIEF EFFECTS CONSIDERED 
EQ. 1, INERITA RELIEF EFFECTS CONSIDERED 
RIGID CONNECTIONS FLAG .......... (RIGID) = 0 
EQ. 0, HINGE CONNECTIONS BETWEEN SOLIDS & SHELLS 
EQ. 1, RIGID CONNECTIONS BETWEEN SOLIDS & SHELLS 


TOTALSYSTEMDATA 


NUMBER OF EQUATIONS .............. (NEQ) = 13 
MAXIMUM HALF BANDWIDTH ............ (MK ) = 1 
MEAN HALF BANDWIDTH.............. (MM) = 8 
NUMBER OF ELEMENTS. .............. (NUME) = 13 
NUMBER OF NODAL POINTS. ............ (NUMNP)= 8 


ORIGINAL NO. OF MATRIX ELEMENTS ........ (NWK) = 34 
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T RESPONSE PRINTOUT (LOADCASE 1) * 
E e e 2 e ICIS e A 3 e e ICIS e ISI e e e e ICICI IC ICICI ICI A I CAI ICA a e e ICA I e ICA a A ae A a fe e a 
xk 


STATUS OF THERMAL EFFECTS FOR THIS LOAD CASE: 


- Thermal flag is off (no thermal effect) 
DISPLACEMENTS 


NODE X-DISPL. Y-DISPL. Z-DISPL. XX-ROT. YY-ROT. ZZ-ROT. 
1 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
2 0.16579 -1.1337 0.00000 0.00000 0.00000 0.00000 
3 0.33157 -0.60429 0.00000 0.00000 0.00000 0.00000 
4 0.38684 -0.33668 0.00000 0.00000 0.00000 0.00000 
5 0.44210 0.00000 0.00000 0.00000 0.00000 0.00000 
6 0.38684  -0.73374 0.00000 0.00000 0.00000 0.00000 
7 0.27631 -0.60429 0.00000 0.00000 0.00000 0.00000 
8 0.16579 -0.33668 0.00000 0.00000 0.00000 0.00000 


MINIMUM/MAXIMUM DISPLACEMENTS 


NODE 1 2 1 1 1 1 
MIN. 0.00000 -1.1337 0.00000 0.00000 0.00000 0.00000 


NODE 5 5 8 8 8 8 
MAX. 0.44210 0.00000 0.00000 0.00000 0.00000 0.00000 


MAXIMUM RESULTANT DISPLACEMENT 


NODE 2 
MAX. 1.1458 


TOTAL STRAIN ENERGY........... = 0.566870E+05 


REACTION FORCE FOR LOAD CASE NO. 1 


NODE CSYS Fx Fy Fz Mx My Mz 
1 0 -.7276E-11 0.7500E-*05 ---------- ---------- ---------- ---------- 
5 0 ---------- 0.2500E+05 ---------- ---------- ---------- ---------- 
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FOR REQUESTED (Global Cartesian Coord. System) 
NODES FX FY FZ MX MY MZ 
Total React. -.7276E-11 0.1000E+06 0.0000E+00 0.0000E+00 0.0000E+00 


0.0000E+00 


SOLUTION TIME LOG IN SEC 


FOR PROBLEM 
TIME FOR INPUT PHASE .............. = 0 
TIME FOR CALCULATION OF STRUCTURE STIFFNESS MATRIX- 0 
TRIANGULARIZATION OF STIFFNESS MATRIX ...... = 0 
TIME FOR LOAD CASE SOLUTIONS .......... = 0 
TIME FOR REACTION/GRID FORCE BALANCE ...... = 0 
TOTAL SOLUTION TIME ..... = 0 


( 0:0: 0) 


STRESS EVALUATION FOR STATIC ANALYSIS 


FKK he e ale ale K K 2K ate ate ate 2K ale ale le ale ale FK 3K ate ate ate FK ale ale K e ale ale 2K FK ate ate ate 2K ale ate K ale ale FK ale 2K ate ate FK FK ale ale 2K ale ale ale ale FK ate ate ate 2K ale ale le le ale ale 2K 2K ate e 


kokod 


£ STRESS PRINTOUT (LOADCASE 1) * 


BAS 8 K K ale ale K 2K ate ate ate ate 2K ale ale K ale ale FK ale ate ate ate FK ale ale FK le ale ale ale ale ate ate ate ate ale ate K ale ale FK ale ale ate ate FK 2K ale ale 2K ale ale ale FK FK ate ate ate ale ale ale le le ale ale 2K 2K FK e 


ek 


STATUS OF THERMAL EFFECTS FOR THIS LOAD CASE: 


- Thermal flag is off (no thermal effect) 


**** NOTE: <STRESS PRINT FLAG IS OFF» **** 


SOLUTION TIME LOG IN SEC FOR STRESS CALCULATIONS 


READING GENERAL INFORMATION AND ELEMENT DATA...= 0 
STRESS CALCULATION AND PRINTOUT ......... = 0 
UPDATING DATABASE ................ = 0 
TOTAL SOLUTION TIME...... 0 

( 0: 0: 0) 


0/09/69/69/9/6/69/69/09/09/9/09/09/09/09/09/09/09/09/09/09/09/09/0/09/0 9/0 0/0/00 0/0/00 00/0/00 %0%0%0 
0/09/69/9/69/69/6/69/09/09/09/09/9/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/0/0/0 


Elemente finite de tip Truss-33 


II.1.2. Aplicatii propuse pentru analiza cu elemente finite de tip TRUSS2D 


Aplicatia II.1.2.1 


Se dà: 


Aa 22500 mm? 


(i =1,2...,11) 


E, 224-10? MPa 


Se cere: 
1 


2. 
3. 
4. 


Aplicatia II.1.2.2 


V 


800 


600 


400 


ODa——F-120KN 


Reactiunile din reazeme 

Fortele axiale din bare 

Tensiunile din bare 

Deplasárile pe orizontalá si verticalá ale punctului A 


Se dà: Aj = ct. 
(i=1,2,...,8) 
A, = 4200mm? 
E, = 24-10? MPa 
v, 20,29 
F=120 KN 


Se cere: 

1. Reactiunile 

2. Tensiunile normale o; din 
bare,i=1,2,....,9 

3 Deplasarea rezultantá din punctul A 
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"E 2 
Aplicatia II.1.2.3 A; = 4700mm 


E, = 2-10% MPa 
vi70,3 
F=250 kN 


Se cere: 
1) Reactiunile 
2) Tensiunile normale o; din 
bare,i=1,2,....,9 
3) Deplasarea rezultantă din 
punctul A 


Aplicatia II.1.2.4 Se da: A; = 20.000mm? 
mM (CERE ce ca 9) 
E, = 21: 10? MPa v;=0,32 
F —400 kN 


Se cere: 

1) Reactiunile 

2) Tensiunile normale o; din 
bare,i=1,2,....,9 

3) Deplasarea rezultantá din punctul A 


Elemente finite de tip Truss-35 


Aplicația II.1.2.5 
Se dá: 
d A; = 1250 mm? 

i ine un 9) 
E - 2.10? MPa 

300 vi=0,3 
F =85 kN 

250 Se cere: 


1) Reactiunile 
| 2) Tensiunile normale o; din 
^ bare,i=l,2,....,9 
3) Deplasarea rezultantá din 
punctul A 


Se dà: 


„= 2 
Aplicația II.1.2.6 dă d Sio 5) 


i A; = 14000 mm” 
(i= 6, 7, 8) 
1 E, = 2-10? MPa 
v;70,3 
(i7 1,2,..., 8) 


F 2185 kN 


Se cere: 

1) Reactiunile 

2) Tensiunile normale o; din 
bare,i=1,2,....,8 

3) Deplasárile pe orizontala si 

verticalá precum si deplasarea 

rezultantá din punctul A 


3 
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Aplicatia II.1.2.7 


c 1250 mm? 
(i= 1,2,...,4) 


i A; = 14000 mm” 
(i — 5, 6, 7, 8) 


E, = 2:10? MPa 
G=1, 2, 3) 
E, =1-10° MPa 
G=4, 5,..., 8) 
Se cere: E 
1) Reactiunile EHE 
2) Tensiunile normale o; din 
bare,i=1,2,....,9 
3) Deplasarea rezultantá din punctul A 
Aplicatia II.1.2.8 
Se dă: 
A; = 12500mm? 
(i= 
152: eoe ,10) 
E = 2.10 MPa 
F =85 kN 
70,3 mm 


Se cere: 

1) Reactiunile 

2) Tensiunile normale o; din 
bare,i=1,2,....,10, ţinând cont de 
imperfectiunea de montaj, A 

3) Tensiunile normale o; din 
bare,i=1,2,....,10, ţinând cont de 
imperfectiunea de montaj, A si 
fortele F aplicate in noduri. 


Elemente finite de tip Truss-37 


Problema: Truss2D-2 

Caracteristici: forte în noduri, variaţie de temperatură, frecvenţe proprii 

Pentru structura plană realizată din bare articulate, şi din acelaşi material, Fig. 
II.1.2, se cere să se calculeze rectiunile, tensiunile din bare şi deplasarea rezultantă a 
punctului A, pentru următoarele cazuri de încărcare: 

a) Încărcarea dată de forțele F; si F3; 

b) Încărcarea dată de o variaţie a temperaturii cu 40 "C faţă de montajul 

initial; 
c) Íncárcare dată atât de forţele F; si F cât si de variaţia de temperatură. 
De asemenea se cere să se calculeze prima frecvență proprie a structurii. 


Fig. II.1.2. Schema geometrică de încărcare si rezemare a structurii plane din bare articulate 


Se dă: Se cere: 
Aj745000 mm” 1) Reactiunile 
GS Decean 10) 2) Tensiunile normale o; din 

E, =1,9-10° MPa bare,i=1,2,....,10 

F,=120 KN 3) Deplasarea rezultantá din punctul A 

F,=150 KN 4) Sá se calculeze aceleasi márimi prezentate 
la punctele 1), 2) si 3), in cazul in care apare 

Densitatea materialului şi o variaţie a temperaturii de A,=40 °C 

p=7,7g/cm=7 7.10% Kg 5) Sá se calculeze frecvențele proprii ale 

° mm: structurii 
Coeficientul de dilatare termicá 
liniará 
PN 
a =11-10° — 


C 
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Rezolvarea problemei cu programul COSMOS/M 
Punctul a) se rezolvă respectând aceleaşi etape prezentate la problema 
Truss2D-1. 


În cele de mai jos se prezintă fişierul de comenzi integral pentru problema 
Truss2D-2. 


VIEW,0,0,1,0 
PLANE,Z,0,1 
PLANE,Z,0,1 
PT,1,0,0,0 
PT,2,700,0,0 
PT,3,2100,0,0 
PT,4,3180,0,0 
SCALE,0 
PT,5,700,1200,0 
PT,6,1500,1200,0 
PT,7,2600,1200,0 
CRPLINE, 1,1,2,3,4,7,6,5,1 
CRLINE,8,5,2 

CRLINE,9,2,6 

CRLINE,10,6,3 

CRLINE,11,3,7 
EGROUP,1,TRUSS2D,0,0,0,0,0,0,0,0 
MPROP,1,DENS,7.7E-6 
MPROP,1,ALPX,11E-5 
MPROP,1,EX,1.9ES 
RCONST,1,1,1,2,45000,0 
M_CR,1,11,1,2,1,1 

NMERGE, 1,22,1,0.0001,0,0,0 
NCOMPRESS, 1,12 
NCOMPRESS, 1,7 
DND,1,UX,0,4,3,UY, 
DND,5,UX,0,5,1,UY, 
DND,7,UY,0,7,1, 
FND,3,FY,-120000,3,1 
FND,2,FY,-240000,2,1 
FND,7,FX,150000,7 


Se activeaza cazul 2 de incárcare, care va 
reprezenta in continuare încărcarea data de 
variatia de temperatura. 


Observatie: 
Cazul 1 de încărcare reprezintă cazul de 
încărcare implicit la deschiderea 
aplicaţiei, deci în cazul nostru reprezintă 
încărcarea dată de forţele aplicate în 
noduri, în conformitate cu secvențele de 
lucru anterioare. 


În mod asemănător în continuare se prezintă celelalte comenzi 


NTCR,1,40,11,1 


Elemente finite de tip Truss-39 


Reprezintă comanda de declarare a încărcării data de 
variatia de temperatura cu 40 °C. 

lata care este calea oferita de meniurile de tip ,,pull down" 
pentru lansarea acestei comenzi: 


d ossi | Control Display Analysis Results Windows Help 


Structural LEI 


Defne by ods 


Fluid Flow ^| Nodal Heat >| Define by Points 
E-Magnetic > Element Heat > 
Load Options >! Heat Flux > Define by Surfaces 
Function Curve ^| Convection ^ | Define by Volumes 
Radiation > Define by Contours 
Hydraulic Flow > Define by Regions 
——————— 
| Delete by Nodes 
Delete by Points 
Delete by Curves 
Delete by Surfaces 
Delete by Volumes 
Delete by Contours 
Delete by Regions 


Plot 
List 


Ca efect se lanseazá fereastra: 


[J NTCR 
Beginning Curve n — | | 
Value [40 — 
Ending Curve m — 
Increment no 


Aceastá fereastrá confine ín ordine: eticheta curbei de 
început a modelului de curbe,1; temperatura aplicată, 

40 "C; eticheta curbei de sfársit a modelului, 11; pasul de 
creştere între etichetele modelului de curbe pe care se 
aplică temperatura. 
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A STATIC,T,0,0,1E-006,1E--010,0,0,0,0,0,0,0,0,0 


Comanda de definire a parametrilor de analiză statică. lata 
care sunt meniurile de tip ,,pull down" pentru lansarea 
acestei comenzi. 


Results ‘Windows Help 


Festart 

Fenumbe- 
Feac:ion 

Data Che-k 

Fun Check 

List Analysis Opticn 


Output Options 


Static 
Frequency/Brcklirg 
Fost Dynamic 
Nonlinear 


Activat> Load Case 
List Lozd Case 
Adaptive Msthcd 
P-Orde- Labels 


Optimize/SensiEiv ty 
Fatique 
Heat_Transfer 
Fluid_Mechancs 
Electro_Magnetic 
Hi-Freq_EMacnetiz 


Stetic Analysis Options 
FFE Stazic Options 
Asvmmstric Load Options 
Stress analysis Options 
PCS Op:ion 


Activate Stress Calc 
Define submodel 
Run Static Analysis 
Run Stress Ana ysis 


Ca efect se lansează fereastra: \<tru-ture 


CJ A STATIC 


Loading flag = 

Inplane effect flag lo: No NS | 
Soft spring flag J 0: No | 

Soft spring value E e-006 

Bcs stiffness value [1 e+010 
Save stiff matrix flag [ 0: No | 
Form stiff matrix flag [ 0: Form | 
Lipdate coordinates flag J 0: No | 
End moment flag for Shell element | 0: Off | 
Grid force balance | 0: Off | 
Inertia relief | 0: Off 
Rigid connections fo: Hinge | 
Solver option [a: Sparse | 
Solution accuracy fo: Off | 


41414 
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A FREQUENCY,1,5,100,0,0,0,0,1 E-005,0,1E-006,0,0,0,0,0,0,0,0,0,0.0 


Comanda de lansare în execuţie a analizei frecvenţelor proprii. Iată 
care sunt meniurile de tip ,, pull down" pentru lansarea acestei comenzi: 


Results Windows Help 


Restart 

Renumber 
Reaction 

Data Check 

Run Check 

List Analysis Option 


Output _Options 
Static 


Frequency Buckling Frequency Options 


Post_Dynamic Buckling Options 
Nonlinear FFE Frequency Option 
Optimizej Sensitivity 
Fatigue 

Heat Transfer 
Fluid_Mechanics 


= Magri 
Ca efect se lansează fereastra: 


CJ A_FREQUENCY 


Number of frequencies [ul 
Method |S: Subspace iteration = 
Maximum number of iterations E oo 

Sturm sequence flag Jo: No bal | 
Eigenvalue shift flag [a: No eigenvalue shift Sol | 

Eigenvalue shift value jo 
Inplane effect flag [ 0: No x | 

Tolerance E e-005 
Soft spring flag Jo: No = | 
Soft spring value E e-006 

Harmonic number jo 

Damping constant [a 
Mass flag [o Lumped mass | 
Modal acceleration flag J 0: No | 
Mass participation factor [ 0: No | 
Save stiffness matrix flag |0: No 
Frequency for nonlinear analysis flag | 0: No j 
Form stiff matrix flag [ O: Form | 
Spin softening flag [ 0: No j 
Rigid connections J 0: Hinge | 
Solver option [o: Spase > = >] 


vv 


Run Frequency 
Run Buckling 


iv" v Y" Y V 4 


CHIE DIE DIE DIE DIE DIE) 


414 
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LCCOMB,51,1,1,2, 


Comanda de lansare in executie a cazului 51 de incarcare 
definit ca fiind obtinut prin combinarea cazurilor 1 si 2 de 
încărcare. 

Iată care este meniul de tip „pull down" pentru lansarea acestei 


comenzi: 
m Windows Help 


Combine Load Case 


Average Modal Stress 
Available Results 
Read Post-Dyn Response 


Setup 
Plot 

List 
Extremes 


CI LCCOMB 


New Load case number 1 MN ^ 

Load case number for term 1 no 3 
Load case factor fortem i | — 
Load case number for term 2 [20 
Load case factor for term 2 ir © 
Load case number for term 3 FO 
Loadcasefactorfortem3/1 = 
Load case number for term 4 
E | AUODJAALLC^ABAOàAE-N 


Notă: 

Cazurile de încărcare cuprinse între limitele 1-50 sunt cazuri 
fundamentale. Cazurile 51-100 sunt rezervate pentru cazuri de 
încărcare obţinute prin combinarea cazurilor fundamentale. 

Pentru un caz de încărcare combinat pot fi combinate maximum 20 
cazuri de încărcare fundamentală, multiplicate însă fiecare cu un 
anumit factor, (care la modul implicit are valoarea 1). 


Elemente finite de tip Truss-43 


Aplicatii: F 


E, 219-10? MPa 
F-20 KN 


Densitatea materialului 


p-2,9g/cm^- 2,9-10% L 


Bae a Tak E mm 
1) Reactiunile : A Coeficientul de dilatare termică 
2) Tensiunile normale o; din bare,i=1,2,....,12 liniară 
3) Deplasarea rezultantă din punctul A d 
4) Sá se calculeze aceleasi márimi prezentate la a —13-10- 9c 
punctele 1), 2) şi 3), în cazul in care apare şi o 
variaţie a temperaturii de A,=25 °C 
5) Sá se calculeze frecvențele proprii ale structurii 
Se dă: 
= 2, - 
Aplicatia II.1.2.10 A=1500mm‘*;  (1-21,2,........ ,8) 
E, 219-10? MPa ; F=2 KN 
22777 Densitatea materialului 
p-7,7g/cm'—7,7 19-5 „ES. 
mm? 
Coeficientul de dilatare termicá 


liniară a 213-107? d 
^C 


1) Reactiunile 

2) Tensiunile normale c; din 
bare,i=1,2,....,8 

3) Sa se calculeze aceleasi marimi 
prezentate la punctele 1), 2) , in cazul 
în care apare şi o variație a 
temperaturii de A,=50 °C 

4) Sa se calculeze frecventele proprii 
ale structurii. 


46 - Elemente finite de tip bará, II 


II.2. Elemente finite de tip bară cu un singur grad de libertate pe nod in 
spatiu (elemente finite de tip TRUSS3D) 


Problema: Truss3D-1 

Caracteristici de încărcare: defecţiuni de montaj 

Să se determine tensiunile care se produc în barele structurii spaţiale, Fig. 
IL.2.1, la montajul forţat al acesteia. A, = 0,5 mm, Ay = 0,6 mm, E = 2.1-10? MPa, aria 
secţiunii transversale a barelor A = 14000 mm”. 


Fig. 11.2.1. Schema geometrică de încărcare şi rezemare a structurii spatiale din bare articulate 


Rezolvare 

1. Definirea geometriei 

2. structurii (puncte şi linii) 
PT, 1,0,0,0 
PT, 2, 0, 1200, 0 
PT, 3, 0, 1200, 1500 
PT, 4, 0, 0, 1500 
PT, 5, 1850, 600, 750 
CRLINE, 1, 1,5 
CRLINE, 2, 2, 5 
CRLINE, 3, 3, 5 
CRLINE, 4, 4, 5 
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I.2, Principalele tipuri de elemente finite-47 


2. Alegerea tipului de element finit, stabilirea proprietátilor de material si a 
constantelor reale 
EGROUP, 1, TRUSSS3D, 0, 0, 0, 0, 0, 0, 0,0 
MPROP, 1, EX, 2.1E5 
RCONST, 1, 1, 1, 1, 14000 
3. Discretizarea curbelor, legarea şi renumerotarea nodurilor 
M_CR, 1, 4, 1, 2, 1, 1 
NMERGE, 1, 8, 1, 0.0001, 0, 1, 0 
Total number of nodes merged = 3 
NCOMPRESS, 1, 7 
Number of Nodes compressed = 2 
4. Definirea conditiilor la limita 
DND, 1, UX, 0, 1, 1, UY, UZ, 
DND, 3, AU, 0, 5, 1, 
DND, 2, UX, 0.5, 2, 1, 
DND, 2, UY, -0.6, 2, 1, 
DND, 2, UZ, 0, 2, 1, 
5. Verificarea datelor de intrare si lansarea in executie 
R_CHECK, STATIC 
R_STATIC 
6. Vizualizarea rezultatelor 
ACTSTR, 1, SX, 2, 1, 0, 0 
STRPLOT, 0, 1, 4, 1, 1, 308.333 


Lin STRESS Le=1 


S1gna X 
62.106 


57.756 
53.406 
43.056 
44,707 
48.357 
36.007 
31.657 


27.397 


Fig. 11.2.2. Tensiunile normale o, din barele structurii reprezentate pentru starea 
deformată 


48 - Elemente finite de tip bará, II 


DISLIST, 1, 2, 1, 5, 1, 0 
Ca efect apare urmatoarea lista care cuprinde fortele de reactiune: 


[_] DISLIST,1,2,1,5,1,0 ES : PIE: 


Node RFX 


1 —3.393e+005 -1.100e+005 -1.376e+005  3.823e+005 
2 2.222e+006 -2.805e+005  0.000e+000 2.240e+006 
3 -7.717e«005 2.503e+005 -3.128e+005 8.695e+005 
4 -7.717e«005 2.503e+005 3.128e+005 8.695e+005 
5 —3.393e+005 -1.100e+005 1.376e+005 3.823e+005 
Sum : 6.250e-002 1.563e-002 0.000e+000 6.442e-002 


Aplicatia II.2.1 


Se dá: 
Y | A=Imm E = 2.109 MPa 

v=0,3 
A172A573A4374A4730.000mm? 
A=Ilmm 
Ay-0,7mm 


Se cere: 


10) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul in care acţionează numai 
forta F 

11) Reactiunile, fortele axiale si tensiunile normale din bare in cazul in care structura este 
solicitată numai de montajul forțat. 

12) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul in care acționează forța F si 
se realizează în acelaşi timp şi montajul forțat; 


I.2, Principalele tipuri de elemente finite-49 


Aplicatia II.2.2 
Se dă: 
E 218-105 MPa E 2188-10? MPa 
A,=A,=A=28.000mm7 A47A3728.000mm? 
A47A471,5A A37A471,5A4 
A,=0,8mm A,=0,2mm 
Ay-0,5mm Ay=0,3mm 


A=0,4mm 

Variatia de temperatura 
A=-30 °C 

Coeficientul de dilatare 
termica a 1.495 1. 


1400 


Se cere: 


7) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul în care structura este 
solicitată numai de montajul forțat. 

8) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul în care structura este 
supusă la variaţia de temperatură fără montaj forțat. 

9) Reactiunile, forţele axiale şi tensiunile normale din bare în cazul în care structura este 
supusă la variaţia de temperatură şi montaj forțat. 


Aplicația II.2.3 Se dă: 
E =1,8-10%MPa 

1900=înălţime A,7A5-28.000mm? 
A37A471 SA 1 
A,=0,2mm 
Ay-0,3mm 
A,=0,4mm 
Variatia de temperatura 
A=-30 °C 
Coeficientul de dilatare 
termica , 44.495 T. 


y=0,3 mm 


50 - Elemente finite de tip bară, II 


Z 
Se cere: 

1) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul în care structura este 
solicitată numai de montajul forțat. 

2) Reactiunile, forţele axiale şi tensiunile normale din bare în cazul în care structura este 
supusă la variaţia de temperatură fără montaj forțat. 

3) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul în care structura este 
supusă la variaţia de temperatură şi montaj forțat. 

Aplicația II.2.4 


Se dă: 

E 22-100 MPa 
A17A5718.000mm? 
A3—2,5A, 

Pentru punctul C: 
A,=0,3 mm 

Variatia de temperatura 
A=50 °C 

Coeficientul de dilatare 
termică 4.15 19-5 L 


Densitatea materialului 
p=7,7g/cm°= 


K. 
Sq 
mm 


Se cere: 

1) Reactiunile, fortele axiale si tensiunile normale din bare in cazul in care structura este 
solicitatá numai de montajul fortat. 

2) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul în care structura este 
supusă numai la variaţia de temperatură. 

3) Reactiunile, forţele axiale şi tensiunile normale din bare în cazul în care structura este 
solicitată numai de forțele indicate. 

4) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul în care structura este 
solicitată prin montaj forțat şi variaţia de temperatură. 

5) Reactiunile, forţele axiale şi tensiunile normale din bare în cazul în care structura este 
solicitată de forțele indicate, variaţia de temperatură şi montajul forțat. 

6) Frecventele proprii înainte şi după montajul forțat. 


I.2, Principalele tipuri de elemente finite-51 


Aplicatia II.2.5 


Se dă: 

E -2.10 MPa 
A17A5718.000mm? 
A3==1,5A, 

Pentru punctul C: 
Ay=2,4 mm 


Variatia de temperatura 
A=30 °C 

Coeficientul de dilatare 
termică 9.15 19-5 as 


Se cere: 


1) Reactiunile, fortele axiale si tensiunile normale din bare in cazul in care structura este 
supusă numai la variația de temperatură. 

2) Reactiunile, forţele axiale şi tensiunile normale din bare în cazul în care structura este 
solicitată numai de forţa F,. 

3) Valoarea forţei F, pentru anihilarea distanţei Ay=2,4 mm faţă de placa rigidă. 

4) Reactiunile, forțele axiale şi tensiunile normale din bare în cazul in care structura este 
solicitată de forțele indicate, variaţia de temperatură şi montajul forțat. 


II.2. Elemente finite de tip bará cu trei grade de libertate pe nod (elemente finite de tip 
BEAM2D) 


Pentru structura cu noduri rigide prezentatá in figurá sá se determine solictárile produse. E — 
2.1-10? MPa, v = 0.3, F = IkN, secțiunea transversală este de formă dreptunghiulară b = 15 mm, 
h = 25 mm. 


Realizarea modelului (preprocesarea): 


1. Pregătirea pentru construirea modelului 
Display>View_Parameter>View 
VIEW, 0, 0, 1, 0 
2. Generarea punctelor (se definesc prin coordonatele lor) 
Geometry>Points>Define 
PT, 1, 0, 0, O 
PT, 2, 250, 0, 0 
PT, 3, 100, 40, 0 
Display>Display_Option>Scale 
SCALE, O 
3. Generarea curbelor 
Geometry>Curves>Line with 2 Pts 
CRLINE, 1, 1, 2 
CRLINE, 2, 2, 3 
CRLINE, 3, 3, 1 
4. Alegerea tipului de elemente finite 
In biblioteca de elemente a programului COSMOS elementul este denumit BEAM2D (Thermal 
and Linear Structural 2D Elastic Beam Element). Este un element finit uniaxial cu douá noduri si 
- trei grade de libertate pe fiecare nod (ux, Uy, rz). 


x 
n group |1 E 


Element Name [is] AM2D: 2D elastic beam element 


Dawe] e x 
OP1:Unused option EO 
OP2:Unused option [a 
OP3:Unused option [a 
DP4:Integration Type [rGass oo x] 
DP5:Material Type [t LinearElastic x] 


PropSets>Element Group OPE:Displacement Formulation [o: Small >] 
EGROUP, 1, BEAM2D, 0, O, O, 0, 0,0,0,0 
OP?:Unused option [a 


OP8:Unused option jo 
Help | Cancel | 
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S. Stabilirea proprietátilor de material 

Pentru analiza structuralá proprietátile de material ce trebuiesc precizate sunt modulul de 
elasticitate E şi coeficientul lui Poisson v. 

PropSets>Material Property 
MPROP, 1, EX, 2.1E5, NUXY, 0.3 

6. Stabilirea proprietăţilor geometrice ale sectiunilor 

Din enunţul aplicaţiei reiese că structura analizată este formată din bare având aceeaşi secțiune 
transversală; aceasta se poate defini alegând din biblioteca programului una din formele 


predefinite, după cum urmează: 
PropSets>Beam Section 
BMSECDEF, 1, 1, 1, 1, 7, 25,15,0,0,0,0,0 


x 
Associated Element group fi 
Real Constant set hi 


Continue Help Cancel 
[Continue ] 
[_] BMSECDEF x] 


Section number JA Rectangular Section 


[ Continue ] Help | Cancel 1 


[_] BMSECDEF x| 
Start location of the section parameters i 


Help | Cancel | 


[_] BMSECDEF x| 
No. of beam section constants to be entered 7 


Help | Cancel | 


x 
BS1:Heightofbeam(H}[25 —— 
BS2:Widlhofbeam (BJ[15 i (iti‘éSC~™S 

BS3 : End-release code [node 1) [D ©. 
BS4 : End-release code [node 2) RN 
BS5 : Shear factor in elem. y-axis [Do ©. 
BSE : Temp. diff. in elem. y-axis [D ©. 
BS7 :Perimeter(HSTAR on |] —tti<“<i‘C;CSC*W 


Help | Cancel | 


End release code este o combinatie de 6 cifre corespunzátoare gradelor de libertate ale nodului. 
Dacá valorile sunt 0 (cazul implicit), valoarea corespunzátoare a fortei nu este cunoscutá, si va fi 
calculată de program. În schimb, dacă în poziţia respectivă este plasat 1, atunci forţa este 
cunoscută ca fiind egală cu zero, iar programul nu va calcula forţa în nodul respectiv. 
7. Discretizarea barelor în elemente finite 
Discretizarea barelor 1 — 3, alegând câte 10 elemente pe fiecare bară: 
Meshing>Parametric_Mesh>Curves 
M_CR, 1, 3, 1, 2, 10, 1 
8. Legarea nodurilor suprapuse si renumerotarea lor 
Meshing>Nodes>Merge 
NMERGE, 1, 33, 1, 0.0001, O, 1, O 
Total number of nodes merged = 3 
Edit>Compress>Nodes 
NCOMPRESS, 1, 32 
Number of Nodes compressed = 2 
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9. Definirea condiţiilor la limită si a încărcării 
Condiţii la limită în nodul ce reprezintă articulația plană din A: 
LoadsBC>Structural>Displacement>Define by Nodes 
DND, 1, UX, 0, 1, 1, UY; 
Condiţii la limită în nodul ce reprezintă rezemarea simplă din B: 
LoadsBC>Structural>Displacement>Define by Nodes 
DND, 11, UY, 0, 11, 1, ; 
Aplicarea fortei F pe verticala: 
LoadsBc>Structural>Force>Define by Nodes 
FND, 21, FY, -1000, 21, 1; 


Solutionarea problemei (procesarea): 


10. Setarea optiunilor pentru analiza tensiunilor si verificarea datelor de intrare 
Analysis>Static>Stress Analysis Options 


A STRESS, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0 x 
Analysis>Run Check AISC code check flac [URTHNINNIEENNNININNINE ~ | 
R_CHECK, STATIC Stress print flag [-1:no,D:all,n:highest n) [o 
11. Lansarea în executie Principal stress print flag E Yes >] 
Analysis>Static>Run Static Analysis Strain print fag [0 No E: 
R STATIC Layer number (post dynamic xy-plots] Boo 


Type of failure analysis Jo: None =] 

Shell face [post dynamic xy-plots) [o: Top >] 
ASME code check faa|0 No x 
J-integral calculation flag [ves = m 

Stress printout at Gauss points fo: No - 
Nodal strain flag [No  —— — rJ 


Help | Cancel | 


Prezentarea rezultatelor analizei (postprocesarea): 
12. Vizualizarea rezultatelor cu ajutorul comenzilor din lista Results 
Tensiunile echivalente Von Mises la nivelul elementelor: 
Results>Plot>Stress 
ACTSTR, 1, VON, 2, 1, 0 
Results>Plot>Stress>Contour Plot 
STRPLOT, O, 1, 30, 1, 1, 2052.52 


Lin STRESS Lc-1 


Von Mises 
10.20100 


9.415600 


8.630200 


7.844800 


7.059400 


6.274000 


5.488500 


4.703100 


3.917700 


Tensiunile echivalente Von Mises la nivelul elementelor 


Diagrama de variație a momentului de încovoiere: 
Results>Plot>Beam Diagrams 
SMPLOT, 1, MT, 1, 30, 1 
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Lin Moment (MT) Lc-1 


Diagrama de variaţie a momentului incovoietor Mr 


19 


Diagrame de eforturi pentru structurá planá de bare cu noduri rigide (BEAM2D) 


F;-4 
Eb (| 24000 A F,- 10000 
El 
E "ym 
O | © 
TAS | Eul APS 
a 60 
E,=2.1E5 ; v=.3 
Ri E;-1.2E5 ; v=.33 


4 — ——29 1000 —_» 4 — — 20— p 


1. Generarea geometriei structurii 
VIEW, O, O, 1, 
PT, 1, O, O, O, 
PT, 2, O, 700, O, 
PT, 3, 1000, O, O, 
PT, 4, 3000, O, O, 
PT, 5, 2000, 700, 0, 
PT, 6, 3000, 700, O, 
SCALE, O, 
CRPLINE, 1, 1, 2, 5, 6, 4, 3, 1, 
CRLINE, 7, 3, 5, 
CRLINE, 8, 3, 6, 
CRLINE, 9, 3, 2, 
2. Alegerea tipului de element finit, stabilirea proprietátilor de material si a constantelor 
reale 
EGROUP, 1, BEAMZ2D, O, O, O, 0, O, O, O, O, 
MPROP, 1, EX, 2.1E5, NUXY, 0.3, 
MPROP, 2, EX, 1.2E5, NUXY, 0.33, 
RCONST, 1, 1, 1, 3, 1200, 3.6E5, 60, 
RCONST, 1, 2, 1, 3, 7153, 4.9E6, 100, 
3. Discretizarea structurii, legarea si renumerotarea nodurilor 
ACTSET, MP, 1, 


ACTSET, RC, 1, 

M CR, 1, 6, 5, 2, 5, 1, 

ACTSET, RC, 2, 

ACTSET, MP, 2, 

M_CR, 2, 3, 1, 2, 5, 1, 

M_CR, 7, 9, 1, 2, 5, 1, 

ACTSET, MP, 1, 

M_CR, 4, 5, 1, 2, 5, 1, 

NMERGE, 1, 54, 1, 0.0001, O, 1, O, 


E 


NCOMPRESS, 1, 53, 1, 

4. Definirea condiţiilor la limită 
DND, 1, UX, 0, 38, 37, UY,, 
DND, 6, UX, 0, 6, 1,, 
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FND, 16, FY, -4000, 16, 1, 
FND, 21, FX, -10000, 21, 1, 
5. Setarea opţiunilor pentru analiza tensiunilor, verificarea datelor de intrare şi lansarea în 
execuţie 
A STRESS, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 
R CHECK, STATIC, 
R STATIC, 
6. Vizualizarea rezultatelor 
SMPLOT, 1, VS, 1, 45, 1, 
SMPLOT, 1, MT, 1, 45, 1, 
SMPLOT, 1, FR, 1, 45. 1, 
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II.3. Elemente finite de tip bará cu sase grade de libertate pe nod (elemente finite de tip 
BEAM3D). 


1. Geometria structurii 
VIEW,0,0,1, 
PT,1,0,0,0, 
PT,2,0,120,0, Fj-300 N 
PT,3,0,300,0, 
PT,4,500,300,0, 
PT,5,500,120,0, 
PT,6,500,0,0, 

SCALE,0, 
CRPLINE, 1,1,2,3,4,5,6,6, 
CRLINE,6,2,5, 


F,=1200 N 


2. Grupul de elemente finite, proprietati de 

material, constante reale 
EGROUP, 1, BEAM3D,0,0,0,0,0,0,0,0, 
MPROP,1,EX,2.1E5,NUXY,0.3 
RCONST,1,1,1,5,800,1.067E5,2.6667E4,40,20, 


3. Discretizarea structurii, legarea nodurilor suprapuse si renumerotarea lor 
VIEW, 1,1,1, 
PT,7,0,0,100, 
M CR,1,2,1,3,8,1,7, 
PT,8,0,120,100, 
PT,9,0,300,100, 
PT,10,500,0,100, 
SCALE,O, 
M CR,3,5,1,3,8,1,9, 
M CR,4,5,1,5,8,1,10, 
M CR,6,6,1,3,8,1,8, 
NMERGE, 1,60,1,0.0001,0,1,0, 
NCOMPRESS, 1,60,1, 


4. Conditii la limitá, incárcare 
DND,1,AL,0,44,43, 
FND,26,FZ,1200,26,1, 
FND,9,FX,300,9,1, 


5. Setarea optiunilor pentru analiza tensiunilor, verificarea datelor de intrare si lansarea in 
executie 

A STRESS, 1, 0, 0, 1, 1, O, O, O, 1, O, O, 

R_CHECK, STATIC, 

R_STATIC, 


6. Vizualizarea rezultatelor 
SMPLOT,1,FR,1,48,1 
SMPLOT, 1, VS, 1,48, 1, 
SMPLOT,1,VT,1,48,1, 
SMPLOT,1,MS,1,48,1, 
SMPLOT,1,MT,1,48,1, 
SMPLOT,1,MR,1,48,1, 
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III. Elemente finite pentru stare planá de tensiune / deformatie / probleme axial-simetrice ( 
III.1. Starea planá de tensiune. 


III.1.1. Placa in stare plană de tensiune (elemente finite tip PLANE2D) / varianta I 


| F\=1000 N 


| F5-2000 N 


grosime ct. = 4,5 mm 


"4 


160 


material OL 52 (E=2.1E5; NUXY= 0.3) 


1. Geometria modelului 
VIEW,O0,0, 1,0, 
PT,1,0,0,0, 
PT,2,0,80,0, 
CRLINE,1,1,2, 
SFEXTR,1,1,1,X,400, 
SCALE,O, 
2. Alegerea tipului de element finit şi setarea corectă a stării plane de tensiune 
PropSets>Element Group 
EGROUP,1,PLANE2D,0,1,0,0,0,0,0,0, 


x| 
Element group |1 E 


Element EnSIPLANE2D: 2D 4-8 node plane element 


[ Continue ] Help | Cancel | 
x 
OP1:Solid/Fluid flag [ata] 
[ Continue | Heb | Cancel | 
— [,) EGROUP X| 
DP2:Integration Type | 1: OM6 >] 
OP3:Type E Plane Stress >] 


OP4:Stress/Strain Direction KAMENA 
1: Axisymmetric 


2: Plane Strain 


3 ET load 


x 


OP5: Material Type OMM aS EE 


OPE: Displacement Formulation | 0: Small m = 
OP7:Material Creep Jo: No Creep >] 
OP8:Strain plasticity Jo: Small >] 


Help | Cancel | 
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3. Definirea proprietátilor de material si a constantelor reale 
MPROP, 1,EX,2.1E5,NUXY,0.3, 
RCONST, 1,1,1,1,4.5, 
4. Discretizarea plácii 
M_SF,1,1,1,4,2,10,1,1, 
5. Conditii la limita si incárcarea plácii 
DCR,1,AL,0,1,1, 
DND,31,UX,0,31,1,UY, 
DND,33,UX,0,33,1, 
FND,9,FY,-1000,9,1, 
FND,18,FY,-2000,18,1, 
FND,32,FX,-2000, 32,1, 
6. Verificarea datelor de intrare si lansarea in executie 
R_CHECK, STATIC, 
R_STATIC, 
7. Vizualizarea rezultatelor 
ACTSTR,1, VON,1,1,0, 
STRPLOT,0,1,20,1,1,1361.39, 
STRLIST,1,1,1,0,1,1,33,1,0, 


Lin STRESS Levi 


Top 
AU XY 


0008-000 —5.638e+000 
.000e+000 —5.299e+000 
.000e+000 —4.960e+000 
.000e+000 -5.170e+000 
.000e+000 —5.299e+000 
.000e+000 —5. 428e+000 
.000e+000 —3.877e+000 
.000e+000 —3.910e+000 
.000e+000 —3.943e+000 
.000e+000 —2.889e+000 
.000e+000 —2.521e+000 
.000e+000 —2.153e+000 
. 000e+000 —2.082e+000 
.000e+000 —2.521e+000 
. 000e+000 —2.961e+000 
.000e+000 2.451e-001 
.000e+000 2.565e-001 


Node 


547e+001 300e+000 
—2.573e+000 -6.568e-001 
2.186e+001 2.987e+000 
-1.817e«001 1.022e+000 
.185e+000 3.408e-001 
1.378e+001 —3.399e-001 

—5.072e+000 —1.724e+000 

-1.698e«000 —3.520e+000 

-2.880e-001 —5.315e+000 
10 5.744e+000 2.523e-001 
11 -2.705e+000 4.647e-001 
12 -9.090e+000 6.771e-001 
13 1.302e«001 -5.500e-001 
14 —2.491e+000 6.726e-001 
15 -1.679e+001 1.895e+000 
16 1.493e+001 —2.518e+000 
17 -1.206e«000 —7.013e+000 


AD CO = S C CO PO IE 
l 
N 


O O O O ODO OOOO OOO O O O OO O O O OO O O O OOO OOOO 


18 —2.128e+001 -1.151e+001 .000e+000 2.680e-001 
19 1.124e+001 -?7.880e-001 .000e+000 2.692e+000 
20 -2.269e+000 6.124e-001 .000e+000 3.034e+000 
21 -1.546e+001 2.013e+000 .000e+000 3.377e+000 
22 2.259e+000 4.737e-001 .000e+000 3.385e+000 
23 —2.246e+000 —1.330e-002 .000e+000 3.034e+000 
24 -6.523e+000 -5.003e-001 .000e+000  2.684e+000 
25 —6.245e+000 -6.551e-001 „000e+000 2.459e+000 
26 —2.920e+000 -1.400e-001 .000e+000 3.034e+000 
27 3.329e-000 3.752e-001 .000e+000 3.609e+000 
28 -6.319e+000  7.244e4000 .000e+000 6.509e-001 
29 -1.180e+001 5.705e+000 .000e+000 3.034e+000 
30 2.115e-001 4.166e+000 .000e+000  5.418e«000 
31 —3.526e+000 —2.193e+001 .000e+000 -3.023e-001 
32 -1.895e+001 -1.696e+001 .000e+000 3.034e+000 
33 3.266e+001 —1.200e+001 .000e+000 6.371e+000 


Distributia tensiunile echivalente Von 
Mises in placa 


CcOcococcooooococococoococcocococoococococcooooooococcoooccorEl 


.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 


Starea plană de tensiuni din placă (tensiuni în nodurile elementelor finite) 


.- ini x! 


TAU YZ 


CcOcocccooooococococococococococococoooooooooococococo 


.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
„000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
„000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
.000e+000 
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Placá in stare planá de tensiune (PLANE2D) / varianta II 


ai 


25 


III.1.2. Placă cu concentrator de tip gaură circulară în stare plane de tensiune (elemente 
finite de tip TRIANG) 


gros.14 mm 


1. Definirea geometriei plăcii 
PLANE,Z,0,1, 
VIEW,0,0,1,0, Fechiv.— 
PT,1,0,0,0, =15000N 
PT,2,150,0,0, 
SCALE, 0, 
PT,3,150,100,0, 
PT,4,0,100,0, 
PT,5,100,50,0, 
CRPCIRCLE, 1,5,3,30,360,4, 
CRPLINE,5,1,2,3,4,1, 


Pregătirea geometriei pentru discretizarea automată (definirea contururilor şi a regiunii) 
CT,1,0,20,4,5,6,7,8,0, 
CT,2,0,15,4,1,2,3,4,0, 
RG,1,2,1,2,0, 
2. Alegerea tipului de elemente finite, definirea proprietatior de material si a constantelor 
reale 
EGROUP, 1, TRIANG, 0,1,0,0,0,0,0,0, 
MPROP,1,EX,2.1E5,NUXY,0.3, 
RCONST,1,1,1,1,14, 
3. Discretizarea automată a regiunii ce reprezintă placa cu gaură circulară 
MA_RG,1,1,1,3,1,0, 
4. Definirea condiţiilor la limită şi a încărcării 
DCR,6,AL,0,6,1, 
FCR,8,FX,-2500,8,1, 
5. Verificarea datelor de intrare si lansarea în execuţie 
R_CHECK,STATIC, 
R_STATIC, 
6. Vizualizarea rezultatelor 
ACTSTR,1,SX,1,1,0,0, 
STRPLOT,0,1,94,1,0, 


Qna X 
32.9100 


28.8740 
24.8370 
20.8010 
16.7650 
12.7280 
8.69190 
4.65550 


9.61913 


Tensiunile normale SX cu evidentierea fenomenului de concentrare a tensiunilor 
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III.2. Starea planá de deformatie. Arc bimetalic solicitat termic (elemente finite de tip 


PLANE2D) 


Arcul bimetalic din figurá este realizat din oţel si aluminiu: Ey = 2.1:10? MPa, vu = 


tensiune la o variaţie a temperaturii AT = 70 °C. 


Material 2 - Aluminiu 


Material 1 - Otel | 


Lungimea 80mm 


1. Definirea geometriei 
VIEW,0,0,1,0 
PT,1,0,0,0 
PT,2,0,0.5,0 
PT,3,0,1,0 
SCALE,O 
CRLINE,1,1,2 
CRLINE,2,2,3 
SFEXTR,1,2,1,X,80 
SCALE,O 


0.3, Qol = 
0.75:10$ PC, Ea; = 7-10* MPa, vai = 0.33, aa; = 12-10 FC. Se studiază starea de deformatie si 


2. Alegerea tipului de elemente finite (setarea corectá a stárii plane de deformatie), 


definirea proprietátior celor douá materiale. 
EGROUP,1,PLANE2D,0,1,2,0,0,0,0,0 


xi 
Element group |1 f 


Element En-SIPLANE2D: 2D 4-8 node plane element 


Ic Help | Cancel | 
E 


= Solid/Fluid flag [iasa 
[ Continue | Help | Cancel " 
x 
DP2:Integration Type fi: OME >] 
OP3:Type [o: Plane Stress >] 
OP4:Stress/Strain Direction 
1: Axisymmetric 


2- Plane Strain 
3: Asisym/asymmetric load 


x 
OPS5:Material Type | ER a a eee 

OP6:Displacement Formulation Er 

OP?:Material Creep [0 Naceep O >] 

OP8:Strain plasticity [o sma OO =] 
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MPROP,1,EX,2.1E5,NUXY,0.3,ALPX,0.75E-6 
MPROP,2,EX,7E4,NUXY,0.33,ALPX,12E-6 
3. Discretizarea celor doua suprafete atribuind proprietati diferite de material. Legarea 
nodurilor suprapuse si renumerotarea lor. 
ACTSET,MP,1 
M_SF,1,1,1,4,8,500,1,1 
ACTSET,MP,2 
M_SF,2,2,1,4,8,500,1,1 
NMERGE, 1,9018,1,0.0001,0,1,0 
NCOMPRESS, 1,9018 
4. Conditii la limitá si definirea incárcárii termice (temperaturá aplicatá in toate nodurile celor 
douá suprafete). 
DCR,1,AL,0,2,1 
LoadsBC»Thermal»Temperature» Define by Surfaces 


NTSFLS021 x 
Beginning Surface Moo 
Vadue[5üD — 
Ending Surface [2 VS 
Increment INESSE 


Help | Cancel | 


5. Setarea opțiunilor analizei pentru analiza statică liniară, verificarea datelor de intrare, 


lansarea în execuţie. 
Analysis>Static>Static Analysis Options 
A_STATIC,T,0,0,1E-006,1E+010,0,0,0,0,0,0,0,0,0,0,0,0 


R_CHECK, STATIC 
ee x 
Loading flag |N: None >] 


Inplane effect flag |N: None 
C: Centrifugal 
Soft spring flag | G: Gravit 
ili: Thermal 
Soft spring value CG: Cent. and Grav. 
Bes stiffness value | CT: Cent. and Therm. 


GT: Grav. and Therm. 
Save stiff matrix flag | CGT: All Three 


Form stiff matrix flag | 0: Form S 


Update coordinates flag [o: No -| 
End moment flag for Shell element | 0: Off bad | 
Grid force balance [rot O] 
Inertia relief | 0: Off z 
Rigid connections | 0: Hinge =] 
Solver option | 0: Sparse m | 
Solution accuracy [rot — e| 
Change to 2nd order | 0: Off >] 
Initial contact clearance | 0: no action -| 
Large displacement [rot | 


Help | Cancel | 


6. Vizualizarea rezultatelor — deplasările UY ale arcului bimetalic si tensiunile normale SX. 
ACTDIS,1,UY,0 
DISPLOT,0,1,8000,1,1,2.3984 
ACTSTR,1,SX,2,1,0,0 
STRPLOT,0,1,800,1,0 
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Lin DISP Lc-1 


Disp Y 


9.99928491 


-0.4167000 


et 


=i 


-22. 


N 


-2. 


Deplasárile verticale UY ale arcului bimetalic 


Lin STRESS Lc-1 


8336800 


„2307009 


„5676009 


0846000 


-5016000 


9186000 


„3355000 


Distribuţia tensiunilor normale SX în cele două materiale 


29 


III.3. Probleme axial-simetrice (element finit de tip axial-simetric PLANE2D). 


Starea de tensiune si deformatie într-un recipient sferic supus la presiune interioară: 


| IModelul 
| real 


| Modelul de calcul 


y Elemente 
finite de tip 
axial - metric 


= 
fff 


Definirea geometriei: 
VIEW,0,0,1, 
PT,1,0,0,0, 
PT,2,100,0,0, 
PT,3,140,0,0, 
CRLINE,1,2,3, 
SFSWEP,1,1,1,Z,90,1, 
Definirea tipului de element finit, a proprietátilor de material si a constantei reale: 
EGROUP, 1,PLANE2D,0,1,1,0,0,0,0,0, 
RCONST,1,1,1,1,1, 
MPROP,1,EX,2.1E5,NUXY, 0.3, 
Discretizarea modelului geometric: 
M SF,1,1,1,4,3,16,1,1, 
Definirea condițiilor la limită pentru problema axial-simetricá şi a încărcării: 
PCR,3,2.5,3,1,2.5, 
DCR,1,UY,0,1,1, 
DCR,2,UX,0,2,1, 
Verificarea datelor de intrare şi lansarea în execuție: 
R_CHECK,STATIC, 
R_STATIC, 
Vizualizarea rezultatelor (distribuţia tensiunilor Von Mises pe grosimea peretelui): 
ACTSTR,1,VON, 1,1,0, 
STRPLOT,0,1,48,1,0, 
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Lin STRESS Lc-1 


Von Mises 
„431809 


a 


5.033100 


4.634300 


4.235600 


w 


„8369300 


3.438209 


3.039500 


N 


„540809 


N 


„242109 


Distribuţia tensiunilor Von Mises pe grosimea peretelui rezervorului sferic 
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IV. Elemente finite de tip placá si invelis 


IV.1. Elemente finite de tip înveliş cu 3 noduri (Shell3) 


Placá plană solicitată la incovoiere 


Presiune pe toata suprafata lateralá a 
plácii de 2 N/cm? 


Realizarea modelului (preprocesarea): 
1. Generarea geometriei modelului 


Definirea punctelor prin coordonatele lor: 

Display>View_Parameter>View 

VIEW, 0, O, 1 
Geometry>Points>Define 

PT, 1,0,0,0 

PT, 2, 200, 0, 0 

PT, 3, 0, 120, 0 

PT, 4, 200, 120, 0 
Display>Display_Option>Scale 

SCALE, 0 
Geometry>Points>Define 

PT, 5, 120, 50, 0 


Generarea curbelor: 
Geometry>Curves>Polyline with Pts 
CRPLINE, 1, 1, 3, 4,2, 1 
CRPLINE, 5, 6, 7,8,6 


Definirea planului XY, înainte de a genera un cerc în acest plan: 


Geometry>Grid>Plane 
PLANE, Z, 0,1 


x 
Axis Normal to Plane [z - | 
Offset on Axis [o 


Grid Line Style [alistan 


Help | Cancel | 


x| 
Curve Poo = 
Keypoint 1 n. - 
Keypoint 2 |  - 
Keypoint 3 Boo 
Keypoint 4 D. ©. 
Keypoint 5 iS 
Keypoint 6 NENNEN 
Keypoint 7 | 5 —— ý 
Keypoint 8 | 
Keypoint 9 |] 1 | 
Keypoint 10[ — tititi‘is~sS 
Kepit I 
Keypoint 12 |] 1 | |— 
Keypoint 13 EN 
Keypoint 14 FY 
Keypoint 15 fe 
Keypoint 16 fe 
Keypoint 17 
Keypoint 18 ; 1 1 | 1] 
Keypoint 19 | 1 ] - 


Help | Cancel | 
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GeomerysCunves»Cices-Cice in Plane TEETER 
CRPCIRCLE, 8, 5, 4, 30, 360, 4 Cuve[8 — —<—s—s 


Keypoint at center 5 © 

Keypoint on zero degree line Moo 
Radius [30 

Angle of the arc [80 | ‘ ‘ ; ;;*;*;éi‘it: 

Number of segments Bo 


Help | Cancel | 


Definirea celor trei contururi ce delimitează placa (contururile se definesc setând drept parametru 


de discretizare dimensiunea elementelor finite): 
Geometry>Contours>Define 
CT, 1, O, 40, 1,1,0,1 
CT, 2, 0, 35, 1,7,0,1 
CT, 3, 0, 35, 1, 10, 0, 1 


n— A 


Contour — 


Mesh flag | 0: Element size >] 


[ Continue ] Help | Cancel | 
Average element size pooo 
Number of reference boundary curves DO O OoOO | 
Help | Cancel | 
Cuvelf. — 
Use selection set Lc 
Redefinition criterion (if curve shared) [ANI TA SIN 


[ ok ] Help | Cancel | 


Definirea unei regiuni delimitată de cele trei contururi, primul exterior iar celelalte două 


interioare: 
Geometry>Regions>Define 
RG, 1, 3, 1, 2, 3,0 


Li xi 
Region i — | | 
Number of contours qT ©. 
te | tee 
Outer Contour is 
Inner Contour 1 - L5 4 
Inner Contour 2 Iia; 1 1 
Undelingsuface[D — 


2. Alegerea tipului de element finit, stabilirea proprietátilor de material si a constantelor 


reale 
PropSets>Element Group 
EGROUP, 1, SHELL3, 0, 0, 0, 0, 0, 0, 0, 0 
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EGROUP 


EGROUP 


EGROUP 


PropSets>Material Property 

MPROP, 1, EX, 2.1E5, NUXY, 0.3 
PropSets»Real Constant 

RCONST, 1, 1, 1, 1, 4.5 


RCONST 


RCONST 


RCONST 


3. Discretizarea automata a regiunii 
Meshing>Auto_Mesh>Regions 
MA_RG, 1, 1, 1, 3, 1, O 
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4. Definirea conditiilor la limitá si a incárcárii (presiune aplicatá pe toatá suprata plácii) 
LoadsBC?»Structural»Displacement» Define by Curves 
DCR, 1, AL, 0, 1, 1, ; 
LoadsBC»Structural» Pressure» Define by Regions 


PRG, 1, 0.02, 1, 1, 0, 4 x 
Beginning Region IEEE 
Pressure magnitude fo. 02 
Ending Region ff 
Increment fi 
Unused option jo 
Pressure direction [tNomd x] 


Help | Cancel | 


Solutionarea problemei (procesarea): 
5. Verificarea datelor de intrare şi lansarea în execuţie 
Analysis>Run Check 
R_CHECK, STATIC 
Analysis> Static>Run Static Analysis 
R_STATIC 


Prezentarea rezultatelor analizei (postprocesarea): 
6. Vizualizarea rezultatelor cu ajutorul comenzilor din lista Results 
Results>Plot>Displacement/Response/Reaction 
ACTDIS,1,UZ,0 
DISPLOT,0,1,46,1,1,8.543 


Disp Z 
„34030 


N 


b 


„04830 


- 


. 75570 


- 


.46300 


- 


„17040 


e 


„87783 


e 


.58522 


e 


„29261 


e 


-00000 


Deplasările plăcii în direcția axei Z (direcția de aplicare a presiunii) 
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Starea de tensiune intr-un element 
de legătură (coltar) pentru două 
cazuri de încărcare. 

E = 2.1:10 MPa, v = 0.3; 


grosime 5mm 


Realizarea modelului (preprocesarea): 

1. Definirea geometriei modelului 
PT,1,0,0,0 
CREXTR,1,1,1,Z,127 
CREXTR,1,1,1,Y,127 
CREXTR,3,3,1,Z,-127 
SCALE,O 
CRFILLET,4,2,1,12.7,1,0,1E-006 
CRFILLET,5,5,2,12.7,1,0, 1E-006 
SFEXTR,1,5,1,X,127 
SCALE,O 
PT,17,63.5,0,69.8 
PLANE, Y,0,1 
CRPCIRCLE,17,17,7,31.75,360,6 
PT,24,19.05,0,31.75 
CRPCIRCLE,23,24,7,6.35,360,6 
CSANGLE, 3,0,63.5,0,69.8,0,0,0,0 
CRGEN,3,23,28,1,1,0,-90,0 
CT,1,0,6.35,3,6,8,1,0 
CT,2,0,6.35,1,22,0 
CT,3,0,3.17,1,26,0 
CT,4,0,3.17,1,32,0 
CT,5,0,3.17,1,36,0 
CT,6,0,3.17,1,42,0 
RG,1,6,1,2,3,4,5,6,0 
ACTDMESH,RG,1 
RGGEN,1,1,1,1,0,0,127,-139.7 
RGSF,2,4,2,6.35 
RGSF,5,5,1,6.35 
CRLINE, 79,15,13 
CRLINE,80,14,11 
CRLINE,81,10,2 
CRLINE,82,4,6 
CT,16,0,6.35,4,47,16,9,79,0,1 
CT,17,0,6.35,4,80,6,15,9,0,1 
CT,18,0,6.35,4,3,5,2,82,0,1 
CT,19,0,6.35,4,81,1,4,2,0,1 
RG,6,1,16,0 
RG,7,1,17,0 
RG,8,1,18,0 
RG,9,1,19,0 

2. Alegerea tipului de elemente finite, definirea constantelor reale si a proprietatilor de 


material 
EGROUP,1,SHELL3,0,0,0,0,0,0,0,0 
RCONST,1,1,1,1,5 
MPROP,1,EX,2E5,NUXY,0.3 
3. Discretizarea modelului, legarea nodurilor şi renumerotarea lor 
MA_RG,1,9,1,3,1,0 
NMERGE, 1,2387,1,0.0001,0,1,0 
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NCOMPRESS, 1,2387 
SHADE, 1,4,1,0.2 
EPLOT,1,3980,1 
RGREORNT,1,3,1 
RGREORNT,7,7,1 
EPLOT,1,3980,1 

4. Definirea condiţiilor la limită si a celor două cazuri de încărcare 
DRG,1,AL,0,1,1 
PRG,2,0.5,2,1,0,4 
ACTSET,LC,2 
PCR,49,-2.5,51,1,-2.5,4 


Solutionarea problemei (procesarea): 

5. Verificarea datelor de intrare şi lansarea în execuţie 
R_CHECK, STATIC 
R_STATIC 


Prezentarea rezultatelor analizei (postprocesarea): 
6. Vizualizarea rezultatelor cu ajutorul comenzilor din lista Results 
ACTSTR,1,VON,1,1,0 
STRPLOT,0,1,3980,1,1,48.3857 
ACTSTR,2,VON, 1,1,0 
STRPLOT,0,1,3980,1,1,1460.47 


Tensiunile echivalente Von Mises 
pentru cazul 1 de incárcare 


Von Mises 
169.4600 


148.2800 
127.1098 
105.9100 
84.73100 
63.54900 
42.36600 
21.18300 


0.000000 


Von Mises 


s.sssss Tensiunile echivalente Von Mises 
78:9» pentru cazul 2 de încărcare 

6.070200 

5.058500 

4.046800 

3.035100 

2.023400 

1.011700 


0.000000 
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IV.2. Elemente finite de tip invelis cu 4 noduri (Shell4) 


Transportor cu melc 


: 2 
presiune p = 0.5 N/mm 
aplicata de-a lungul melcului pe un sfert din 


suprafata laterala a spirei 


grosime 2.5mm 


pasul = 60 mm 


1. Generarea geometriei modelului 
PT,1,40,0,0, 
PT,2,80,0,0, 
PTGEN,1,1,2,1,2,0,0,15,0,0,90, 
PT,5,0,0,0, 
CRARC,1,1,3,5,41.3824, 
CRARC,2,2,4,5,80. 7001, 
SF2CR,1,1,2,0, 
SFEXTR,1,1,1,Z,60, 


Generarea geometriei melcului prin roto-translatie asupra suprafeţelor 1 şi 2: 
Geometry>Surfaces>Generation>Generate 

SFGEN,11,1,2,1,2,0,0,15,0,0,90, 

SCALE,O, 


xi , 
Generation number nt | 
Beginning Surface NENNEN : 
Ending Surface [2 Vc 
Increment NEN x 


Generation flag Amii 


Help | Cancel | 


x 
X-Displacement [III 
Y-Displacement NN 
Z-Displacement is ©. 

X-Rotation jf © | 
Y-Rotation Doo 
Z-Rotation [uà ^. 


Help | Cancel | 
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Alegerea tipului de elemente finite, definirea constantelor reale si a proprietátilor de 


material 
EGROUP,1,SHELL4,1,0,0,0,0,0,0,0, 
RCONST,1,1,1,1,2.5, 
MPROP,1,EX,2.1E5,NUXY,0.3, 


3. Discretizarea structurii, legarea nodurilor suprapuse si renumerotarea lor 
M_SF,1,24,1,4,2,2,1,1, 
NMERGE, 1,216,1,0.0001,0,1,0, 
Total number of nodes merged=108 
NCOMPRESS, 1,216, 
Number of Nodes compressed-108 


4. Definirea condiţiilor la limită si a încărcării 
DCR,48,ALL,0,54,2, 
PSF,1,0.5,17,8,0.5,0.5,4 


Structura după discretizare şi 
aplicarea condiţiilor la limită 


5. Verificarea datelor de intrare şi lansarea în execuţie 
R CHECK, STATIC, 
R STATIC, 
6. Vizualizarea rezultatelor 
ACTSTR,1, TXY,1,1,0,0, 
STRPLOT,0,1,96,1,1,25.686 


Lin STRESS Lc-1 


Teu XY 
52.8450 


41.5600 
30.2740 
18.9880 
7.79180 
-3.5841 
-14.870 
-26.156 


-37.442 


Tensiunile tangentiale TXY produse în structură 
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V. Elemente finite tridimensionale 
V.1. Elemente finite hexaedrice (elemente finite de tip SOLID) 


Să se determine starea de tensiune în bara curbă volumicá din figură. E = 2.1-10? MPa, v 
= 0.3, F) = 3kN, F - 2 KN. 


F;-3000 N 


F5-2000 N 


E-2.1E5 MPa 
v=0.3 


Realizarea modelului (preprocesarea): 


1. Generarea geometriei modelului 
Geometry>Points>Define 
PT, 1, 0, 0,0 
PT, 2, 120, 0, 0 
PT, 3, 120, 20, 0 
Geometry>Curves>Line with 2 Pts 
CRLINE, 1, 2, 3 
Generarea prin extrudare a suprafetei ce reprezintá sectiunea transversalá dreptunghiulará a 


barei: l 
Geometry>Surfaces>Generation>Extrusion xi 
SFEXTR, 1, 1, 1, X, 20 Beginning Curve E 
Ending Curve E 
Increment E 
Axis symbol [x - | 
Value [20 
Generarea volumului prin rotirea sectiunii transversale Hep | Cancel | 
dupá axa y: 
Geometry» Volumes»Generation»Sweeping x| 
VLSWEEP, 1, 1, 1, Y, -90, 1 
Beginning Surface E 
Ending Surface E 
Increment E 
Axis symbol [v - | 
Angle of the arc [:30 
Number of segments E 


Help | Cancel | 40 


Definirea unui sistem de coordonate cartezian a cárui origine sá coincidá cu punctul de aplicatie 
al forţei F», iar direcția X cu direcţia acestei forte: 
Geometry>Curves>Manipulation>Break (equally) 


iii ea E 
Beginning Curve | 
Ending Curve no 
Increment N] 
Number of curve segments [| V 


Original curve keeping flag [39 BT a 4-15 


Help | Cancel | 


Geometry»Coordinate Systems»3 Points 
CSYS, 3,0, 11,6,5 x 
Coordinate system [D ©. 
Coordinate system type [0 Cartesian x] 
Keypoint at origin nos o 
Keypoint on the X-axis ff © 
Keypoint on the X-Y plane fi o 


Help | Cancel | 


2. Alegerea tipului de elemente finite 
Se alege din biblioteca progamului elementul finit de tip SOLID (3D 8- to 20-Node 


Isoparamatric Solid: 
PropSets>Element Group 
EGROUP, 1, SOLID, 0, 1, 0, O, O, 0, 0, 0 


xi 
Element group E 


Element Name $m N] DA 20-node 3D solid element 


x] Help | Cancel | 
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Tjecrour 0 


OP1:Solid/Fluid flag [eR] 


Help | Cancel | 


x 
DP2:Integration Type 
OP3:Unused option ae 
OP4:Stress/Strain Direction [Global  — rl 
ÜP5:Material Type JO: LinearElastie x] 
OP6:Displacement Formulation [Smi o 
OP?:Material Creep [r NoGep O >] 
OP8:Strain plasticity [rsmd — — rl 


Help | Cancel | 


3. Definirea proprietátilor de material 
Pentru analiza structuralá proprietátile de material ce trebuiesc precizate sunt modulul de 
elasticitate E şi coeficientul lui Poisson v. 
PropSets>Material Property 
MPROP, 1, EX, 2.1E5, NUXY, 0.3 


4. Discretizarea volumului 
Discretizarea volumului se face parametric, utilizatorul alegând numărul de elemente finite pe 


fiecare dintre cele trei directii ce definesc laturile volumului. 
Meshing»Parametri Mesh» Volumes 
M VL, 1, 1, 1, 8, 2, 2,12, 1, 1,1 


xi 
Beginning Yolume Moo 
Ending Volume ho 7 
Increment |1 NENNEN 
Number of nodes per element [B sr 
Number of elements on first curve Poo 
Number of elements on second curve [2 € 
Number of elements on third curve ECE 
Spacing ratio for first curve Moo 
Spacing ratio for second curve Moo 
Spacing ratio for third curve i 


Help | Cancel | 


5. Definirea condiţiilor la limită si a încărcării 
Condiţii la limită pe suprafaţa ce reprezintă încastrarea şi de-a lungul laturii rezemate: 
LoadsBC>Structural>Displacement>Define by Surfaces 


571051: BNET s 


Beginning Surface n | 


Displacement label cae All translations >] 


Value Doo 


Ending Surface f 
Increment 
Additional Displacement labels if any [L2.....L6,) | 


Help | Cancel | 
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Loads»Structural»Displacement» Define by Curves 


DCR, 8, AU, 0, 8, 1,; 
— | 


Beginning Curve BB ©. 


Displacement label Aer All translations ~] 


Value NENNEN 


Ending Curve [8 
Increment fi 
Additional Displacement labels if any (L2.....L6,) | 


Help | Cancel | 


Definirea încărcării reprezentată de cele două forte: F; — a cărei direcție se defineşte în sistemul 
de coordonate global (denumit sistemul de coordonate 0) şi Fz — a cărei direcție se defineşte în 
sistemul de coordonate cartezian denumit 3. 
Control>Activate>Set Entity 
ACTSET, CS, 0 
LoadsBC>Structural>Force>Define by Nodes 
FND, 111, FX, 3000, 111, 1 
Control>Activate>Set Entity 
ACTSET, CS, 3 
LoadsBC>Structural>Force>Define by Nodes 
FND, 63, FX, 2000, 63, 1 
Nodurile in care s-au aplicat cele douá forte (nodul 111 si respectiv 63) au fost indicate cu 
mouse-ul (click stánga). 


Solutionarea problemei (procesarea): 


6. Verificarea datelor de intrare si lansarea in executie 
Analysis>Run Check 
R CHECK, STATIC 
Analysis>Static>Run Static Analysis 
R_STATIC 


Prezentarea rezultatelor analizei (postprocesarea): 
7. Vizualizarea rezultatelor cu ajutorul comenzilor din lista Results 
Results>Plot>Stress 
ACTSTR, 1, VON, 1, 1, 0, 0 
Results>Plot>Stress>Contour Plot 
STRPLOT, 0, 1, 48, 1, 1, 652.033 


Von Mises 
161.2100 


141.3188 
121.4200 
101.5200 
Y 81.62500 
61.72900 
41.83300 
21.93700 


2.041600 


Distribuția tensiunilor echivalente Von Mises reprezentată pe deformata modelului 
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V.2. Elemente finite tetraedrice (elemente finite de tip TETRA10) 


Discretizare automată 3D a unui dispozitiv cu 
deschideri 


1. Definirea geometriei 
PLANE,Z,0,1 
VIEW,0,0, 1,0 
Construirea feței laterale dreptunghiulare (fata din dreapta în figura): 
GRIDON,0,0,5,5,20,20,2 
CRPCORD,1,0,50,0,0,0,0,80,0,0,80,50,0,0,50,0 
Rotunjirea colturilor superioare: 
CRFILLET,5,1,4,8,1,0,1E-006 
CRFILLET,6,3,4,8,1,0,1E-006 
GRIDOFF 
PTDEL,5,8,3 
PTCOMPRESS,1,10 
REPAINT 
Construirea a două cercuri, ce reprezintă urma deschiderilor transversale cilindrice în 
planul acestei feţe laterale: 
PT,7,20,35,0 
PT,8,60,35,0 
CRPCIRCLE, 7, 7,3,10,360,4 
CRPCIRCLE,11,8,5,10,360,4 
Transformarea dreptunghiului cu cercuri într-o regiune (la definirea contururilor 
dimensiunea elementelor s-a definit ca fiind egală cu unitatea, aceasta urmând a fi redefintă la 


construcția poliedrului): 


CT,1,0,1,1,3,0 
CT2014180 4 Mew pun E 
CT,3,0,1,1,14,0 " N N 
RG,1,3,1,2,3,0 
VIEW,1,1,1,0 
ho d LEM T 
Contururi 1, 2, 3 
1 
R1 
Regiunea R1 


Construirea modelului 3D prin rabateri si extrudări: 
RGGEN, 1,1,1,1,0,0,0,60 
SFEXTR,4,14,1,Z,60 
SFEXTR,2,2,1,Z,60 
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Regiunile R1, R2 
Suprafetele superioare 1, 2, 3 
Suprafetele cilindrice 4 -11 
Suprafata inferioará 12 


Constructia deschiderii longitudinale (urma acesteia in planul fetei laterale este un 
dreptunghi): 

VIEW,1,0,0,0 
PLANE,X,0,1 
GRIDON,0,0,5,5,20,20,2 
CRPCORD,43,0,20,10,0,20,50,0,5,50,0,5,10,0,20,10 
GRIDOFF 
VIEW,-3,2,1,0 
REPAINT 
CT,7,0,1,1,46,0 
CT,8,0,1,4,19,41,1,31,0,1 
RG,3,2,8,7,0 
RGGEN,1,3,3,1,0,80,0,0 
SFEXTR,43,46,1,X,80 
REPAINT 
CLS,1 


Regiunile R3, R4 
Suprafeţele deschiderii longitudinale 13 -16 


Generarea entităților geometrice Polyedra (cu redefinirea dimensiunii elementelor finite) 
şi Part pregătitoare pentru discretizarea automată 3D 
Geometry>Polyedra>Define 
PH,1,SF,1,10,0.0001,1 


de o x 
Polyhedron label Moo 
Reference entity type 
Help | Cancel | 
xi 
Surface n. | | 
Average element size ho 9 € 
Tolerance ponn 


Redefine element size on boundary 3E zs 


Help | Cancel | 
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2. Definirea grupului de elemente finite si a proprietátilor de material 
EGROUP,1,TETRA10,0,1,0,0,0,0,0,0 
MPROP,1,EX,2.1E5,NUXY,0.3 
3.Discretizarea automata a părții 
MA_PART,1,1,1,1,1,4,0,0 
4. Stabilirea condiţiilor la limită şi a încărcării: dispozitivul este fixat pe faţa sa inferioară şi 
supus la presiune normală pe suprafeţele cilindrice interioare (p = 5 Mpa). 


DSF.124L0124 x 
PSE AS 125,34 Beginning Surface ao ©. 
Pressure magnitude 56 © 

Ending Surface LEES 

Increment NENNEN 

Pressure at the end of direction 1 56 V 

Pressure at the end of direction 2 bt © 

Pressure direction [tNomd — | 


Help | Cancel | 


5. Verificarea datelor de intrare şi lansarea în execuţie 
R_CHECK,STATIC 
R_STATIC 
6. Vizualizarea rezultatelor 
ACTSTR,1,VON,1,1,0 
STRPLOT,0,1,1269,1,0 


Tensiunile Von Mises produse în 
dispozitiv 
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VI. Analiza câmpurilor termice 
VI.1. Tensiuni termice într-un tub cilindric 


Se analizează distribuția temperaturii şi a tensiunilor termice pe grosimea peretelui unui 
tub cilindric din oţel, temperatura interioară fiind Tint = 80 °C, iar cea exterioară Text = 20 °C. 
R = 20 cm 
r=10cm 
Text = 80 °C 
Tint = 20 °C 


1. Definirea geometriei 
VIEW, O, 0, 1, 0 
PT, 1, 0.1, 0, 0 
PT, 2, 0.1, 0.01, 0 
CRLINE, 1, 1, 2 
SCALE, O 
SFEXTR, 1,1, 1, X, O.1 
2. Definirea grupului de elemente finite (PLANE2D pentru probleme axial-simetrice) si 


alegerea materialului din biblioteca programului 
EGROUP, 1, PLANE2D, 0,1,1,0,0,0, 0,0 
PICK MAT, 1, STEEL, SI 
MPLIST, 1, 1, 1 


[]MPLIST,1,1,1 i 


Label Name Temp/BH_Cr Value 
À 1 n 2.000000e+011 
A 1 NUXY D 3.200000e-001 
A 1 GXY 0 ?.600000e+010 
A 1 ALPX 0 1.200000e-005 
À 1 DENS D ?.800000e+003 
A 1 C 0 5.000000e+002 
A 1 KX n 3.000000e+001 
A 1 MPERM R 0 1.000000e+000 


3. Discretizarea suprafeței 
M_SF, 1, 1, 1, 4, 2, 20,1, 3 


4. Definirea condițiilor la limită corespunzătore xj 
problemei axial-simetrice şi a temperaturii pe Begirning Curve | 


exteriorul şi interiorul cilindrului tubular. 


DND, 1, UY, 0, 63, 1, UZ, RX, RY, RZ Value [353 

ACTSET, TC, O Ending Curve fi 

NTCR, 1, 353, 1, 1 OO —— ————— 
NTCR, 2, 293, 2, 1 

ACTSET, LC, 1 Hep | _ Cancel | 
ACTSET, TC, 1 


57 


D 3 
q 
Å 
Â 
A A 


AA 


Modelul de calcul 


5. Setarea opţiunilor pentru analiza termică şi lansarea în execuţie 
Analysis>HEAT_TRANSFER>Thermal Output Options 
HT_OUTPUT, 1, 1, 1 
Analysis>HEAT_TRANSFER>Thermal Analysis Options 
A THERMAL, S, 0.001, 1, 1, 20, 0, 1, 1E+008, 0, 0, O 
Analysis>HEAT_TRANSFER>Run Thermal Analysis 
R_THERMAL 


6. Vizualizarea rezultatelor (variatia temperaturii pe grosime) 
Results, Plot, Thermal 
ACTTEMP, 1, TEMP 
Results, Plot, Thermal, Contour Plot 
TEMPPLOT, O, 1, 40, 1, 0 


THERMAL Step=1 


Temp 
353.0 


345.5 
338.0 
330.5 
323.0 
315.5 
308.0 
300.5 


293.0 


Distributia temperaturii in cilindrul tubular [°K] 


7. Definirea temperaturii ca incárcare si analiza tensiunilor termice 
LoadsBC>Loads_Options>Read Temp as Load 
TEMPREAD, 1, 1, 1, 1, 1 
Analysis>Static>Static Analysis Options 
A STATIC, T, 0, 0, 1E-006, 1E+010, 0, O, 0, 0, O, O, O, O, O, O, 0, 0 
Analysis>Static>Run Static Analysis 
R STATIC 


8. Afişarea tensiunilor termice radiale (SX) şi circumferentiale (SZ) 
ACTSTR, 1, SX, 1, 1,0 
STRPLOT, O, 1, 40, 1, 0 
ACTSTR, 1, SZ, 1, 1, 0 
STRPLOT, O, 1, 40, 1, 0 
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Lin STRESS Lc=1 


S1gma X 
-1.6481E+006 


-3.6682E+006 
-5.6884E+006 
-7.7086E+006 
-9.7288E+006 


-1.1749E+007 


.37698E*907 
89E*097 


10E+007 


y 
Eu 


Lin STRESS Lc-1 


Tensiuni termice radiale o, (SX-10° MPa) 


Signa Z 
8.14130E+007 


5.49548E*997 
2.84950E+007 
2.03610E+006 
-2.4423E+007 


-5.0882E+007 


—7 F341 E+007 


Tensiuni termice circumferentiale og (SZ:10* MPa) 
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VI.2. Conductie nestationará într-un cilindru de lungime mare 


Un cilindru lung din aluminiu de diametru 5 cm şi cu o temperatură inițială de 200 ?C, 
este expus brusc într-un mediu de convecţie a cărui temperatură este de 70 *C. Se determină 
temperatura cilindrului după un minut de la expunere. 


raza cilindrului r=0.025 m 
temperatură inițială AL To=200 °C 
temperatura mediu T.=70 °C 
conductivitate termica AL K=215 W/m°C 
densitate AL p-2700kg/m* 
cáldurá specificá AL C=936.82 J/kg °C 


coeficient de transfer termic prin convectie 
h=525 w/m? °C 


Pentru cá geometria si conditiile la limita sunt axial-simetrice, se vor utiliza elemente 
finite Plane2D pentru probleme de tipul axial-simetric. 


Definirea geometriei 
VIEW, O, 0, 1, 0 
SFA4CR, 1, 0, O, O, O, 0.00125, 0, 0.025, 0.00125, 0, 0.025, 0, 0 
SCALE, O 
Definirea grupului de elemente finite 
EGROUP, 1, PLANE2D, 0, 1,1, 0, O, 0, 0, 0 
Definirea proprietátilor de material pentru aluminiu 
MPROP, 1, KX, 215 
MPROP, 1, DENS, 2700 
MPROP, 1, C, 936.82 


Discretizarea modelului geometric 
M_SF, 1, 1, 1, 4, 1, 20, 1, 1 
Definirea conditiilor la limita (convectie termica) 
LoadsBc>Thermal>Convection>Define by Curves 


GG e ai et x 
Beginning Curve 2 - 
Convection coefficient |525 
Ambient temperature HER 
Ending Curve e © 
Increment fi 
Time curve for ambient temperature [o ©. 


Help | Cancel | 
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Definirea temperaturii cilindrului din aluminiu şi a numărului de paşi — 60- pentru analiză 


Loads8C>Load_Options>Uniform Temp x 
TUNIF, 200 : pud 
LoadsBC»Load Options»Time Parameters Uniform temperature |20 
TIMES, 0, 60, 1 Hep [| Cancel | 


x 
Starting Time NN 
FinalTimef60 — 8 || 
Time Increment fi ©. 


Help | Cancel | 


Setarea optiunilor pentru analiza termicá, lansarea in executie 
A THERMAL, T, 0.001, 1, 1, 20, O, 1, 1E+008, 0, 0, 0 
HT OUTPUT,1,1,1 
R THERMAL 


Afişarea rezultatelor (temperatura în cilindru după un minut /ultimul pas) 
ACTTEMP, 60, TEMP 
TEMPPLOT, O, 1, 20, 1 


THERMAL Step=60 


Temp 
119.86 


119.68 
119.49 
119.30 
119.12 
118.93 
118.75 
118.56 


118.37 


Temperatura in cilindru dupa un minut de expunere (°C) 
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V.3. Modelare geometrică SolidWorks si elaborare model calcul 3D CosmosWorks 


ANALIZA STĂRII DE TENSIUNE SI DEFORMATIE DINTR-UN 
ELEMENT DE FIXARE 


1. Scopul lucrării. 

Lucrarea are ca scop realizarea unui element de fixare cu ajutorul programului CAD 
SolidWorks şi rularea unei analize prin Metoda Elementului Finit implementată în 
programul COSMOSWorks, în urma aplicării unei incárcári asupra piesei. 


2. Realizarea elementului de fixare. 

Se initializeazá programul CAD SolidWorks. Pentru crearea piesei se lansează 
comanda File — New — Part din bara de meniu superioară. Din fereastra meniu stângă 
se selectează planul de bază Right Plane . (Fig. 1) 


$i SolidWorks 2005 - [Part1] 


$ Fie Edit view Insert Tools COSMOSWorks Window Help Em 
JD SHB 9-c-[RwE 8/EBm-5bizo ||$jqaqseeg-I]oogsiamc|a-eusmesem 
IETTTEFLELINESTIITIIITITIIZETLIE 


e 
(=) Să) 
e EN 
e o 
e o 
2a 3 
& E 
@ a 
= = 
@ i 
a N 
icy * 
$38 E" 
e & 
a e- 
w- A 
e m 
3 


Lă 
[T] (Editing Part [| 7 


Fig. 1. 


Având planul de bază Right Plane selectat, se lansează comanda Insert — Sketch din 
bara de meniu superioară. Din originea sistemului de axe se schiteazá liber, fără 
dimensiuni, conturul închis ca în figura 2, utilizând succesiv comanda Line , aflată în 
bara de meniu dreaptă. Dimensiunile secțiunii T întors ce definesc conturul închis sunt 
următoarele: - talpa - 30 x 10 mm ; - inima - 25 x 10 mm. 
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Fig. 2. 


Dimensiunile se declará selectánd fiecare entitate in parte (linie), succesiv prin 
intermediul comenzii Smart Dimension aflată in bara de meniu dreaptă şi introducând pe 
rând valorile. Secţiunea transversală T întors, dimensionată, este redată în figura 3. 
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Odatá realizatá sectiunea transversalá, se lanseazá comanda Rebuild din bara de 
meniu superioará, pentru definitivare. Din acest moment sectiunea transversalá a piesei 
apare în bara de meniu stângă sub numele de Sketch1 . Selectând entitatea Sketch1, se 
lansează comanda Insert — Boss/Base — Extrude din bara de meniu superioară. Astfel 
se realizează extrudarea secțiunii transversale. Lungimea de extrudare este de 50 mm. 
Entitatea geometricá obtinutá este redatá in figura 4. 


Tie (de Wem inet Tooke COSMOSWed Widow Map aleiz 
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2t 
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Fig. 4. 


În etapa următoare se selectează cele două muchii înguste aflate la nivelul superior al 
inimii secțiunii T întors şi se lansează comanda Insert — Features — Chamfer din bara 
de meniu superioară. Lungimea pe care se efectuează teşirea celor două muchii este de 25 
mm. În urma operaţiunii de teşire piesa dobândeşte forma ilustrată în figura 5. 
Următoarea fază este materializată în rotunjirea capătului superior al piesei, la nivelul 
prismei triunghiulare rezultată din operaţiunea de teşire. Pentru rotunjire se selectează 
muchia îngustă superioară a prismei triunghiulare şi se lansează comanda Insert — 
Features — Fillet/Round . Raza de rotunjire are valoarea de 9 mm. În urma operaţiunii 
de rotunjire, piesa are alura redată în figura 6. 

În continuare se realizează partea cilindrică a piesei. Pentru aceasta se selectează fata 
orientată spre observator a prismei triunghiulare rotunjite şi se lansează comanda Insert 
— Sketch din bara de meniu superioară. Având fata selectată se trasează liber fără 
dimensiuni un cerc utilizând comanda Circle din bara de meniu dreaptă. Prin intermediul 
comenzii Smart Dimension şi având selectat cercul trasat anterior se declară valoarea de 
18 mm ce reprezintă diametrul. Se selectează cercul dimensionat împreună cu sectorul 
circular orientat spre observator ce aparţine prismei triunghiulare rotunjite şi se alege 
opțiunea Concentric , ce apare în meniul deschis în partea stângă. Ca o ultimă operaţiune, 
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se lanseazá comanda Rebuild din bara de meniu superioará, pentru definitivare. Piesa 
împreună cu cercul trasat este redată în figura 7. 


Fig. 6. 
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Fig. 7. 


Din fereastra stângă se selectează entitatea Sketch2 (ce reprezintă cercul trasat 
anterior) şi se lansează comanda Insert — Boss/Base — Extrude din bara de meniu 
superioară. Lungimea de extrudare este de 5 mm. După această operaţiune, piesa 
dobândeşte alura redată în figura 8. 


Fig. 8. 
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Următoarea operaţiune o constituie gáurirea părţii cilindrice a piesei. Într-o primă fază 
se selectează suprafaţa circulară a părții cilindrice a piesei, orientată spre observator şi se 
lansează comanda Insert — Sketch din bara de meniu superioară. Având suprafața 
amintită selectată se trasează liber fără dimensiuni un cerc utilizând comanda Circle din 
bara de meniu dreaptă. Prin intermediul comenzii Smart Dimension şi având selectat 
cercul trasat anterior se declară valoarea de 8 mm ce reprezintă diametrul găurii. Se 
selectează cercul dimensionat împreună cu cercul orientat spre observator ce aparține 
părții cilindrice a piesei şi se alege opţiunea Concentric , ce apare în meniul deschis in 
partea stângă. Ca o ultimă operaţiune, se lansează comanda Rebuild din bara de meniu 
superioară, pentru definitivare. Piesa împreună cu cercul trasat anterior ce serveşte 
operațiunii de găurire este redată în figura 9. 


fs SclidWorks 2005 - [L1 *] lar 
fb fis (de Ve încat Tooke COGOMOSWeds Wedon Hip igi xi 


JD EE S S9 C [SE 8 DR eH Le ISRAI tT- 00083 aaa nu ana: 
"EI PT 


TETY MAANTIET: 


999&3560608?99 


s4 Be 


Randy Distance between corde oies: Sile Fly Oaled | dking Piet (T i ki 
Fig. 9. 


Din fereastra stângă se selectează entitatea Sketch3 (ce reprezintă cercul destinat 
operaţiunii de găurire) şi se lansează comanda Insert — Cut — Extrude din bara de 
meniu superioară şi se alege opţiunea Through All . În final se lansează comanda 
Rebuild din bara de meniu superioară şi astfel generarea geometriei piesei este finalizată, 
alura ei fiind redată în figura 10. 


3. Definirea şi rularea analizei cu elemente finite. 

Din bara de meniu superioară se lansează comanda COSMOS Works — Study... . În 
fereastra ce se deschide se denumeşte studiul - Study name - (se introduce Studiu 1), se 
selectează opțiunea Static - tipul analizei (Analysis type) - deoarece se consideră că 
sarcinile exterioare se aplică static asupra piesei - şi se alege tipul Solid al discretizării 
geometriei piesei (Mesh type). Aceste faze succesive constituie definirea analizei cu 
elemente finite. 
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Fig. 10. 


Ín continuare se alege materialul din care se considerá realizarea piesei. Pentru aceasta 
se lanseazá comanda COSMOSWorks — Apply Material to All... din bara de meniu 
superioará si astfel apare fereastra de declarare a materialului, ce este redatá in figura 11. 


Material 


1023 Carbon Steel Shee 
201 Annealed Stainless 


AISI 1015 Steel, Cold Dr 
AISI 1020 

AISI 1020 Steel, Cold Rc 
AISI 1035 Steel" 

AISI 1045 Steel, cold dr: 


Launch] 


4286 Iron Base Superalloi 
AISI 1010 Steel, hot rolle | 


| 


i 


Linear Elastic Isotropic El 


Elastic modulus 


Poisson's ratio 


Shear modulus 
Mass density 


Tensile strength 


Compressive strength 


Yield strength 
Thermal expansion coe 


. Thermal conductivity 


Specific heat 


EE 


Fig. 11. 


Se aleg optiunile From library files — cosmos materials — 201 Annealed Stainless 
Steel . Materialul ales reprezintá un otel inoxidabil. 
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În următoarea etapă se definesc condiţiile de rezemare şi încărcare a piesei. Se 
consideră că piesa este fixată la nivelul tălpii inferioare şi deci această suprafață nu se 
poate deplasa sau roti în raport cu axele sistemului de cordonate Oxyz. Pentru a aplica 
această condiţie, în primul rând se roteşte piesa într-o poziţie convenabilă astfel încât să 
fie vizibilă suprafața tălpii inferioare a piesei. Se selectează această suprafață şi se 
lansează comanda COSMOSWorks —» Loads/Restraint — Restraints... din bara de 
meniu superioará. Ín fereastra ce se deschide in partea stángá se alege optiunea Fixed . 
Fixarea tálpii inferioare a piesei este redatá in figura 12. 
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Fig. 12. 


În continuare se aplică încărcarea asupra piesei. Se consideră cá încărcarea se 
materializează într-o presiune în valoare de 50 MPa aplicată asupra feţei găurite a părții 
cilindrice a piesei. Pentru aplicarea încărcării, se selectează fața amintită şi se lansează 
comanda COSMOSWorks — Loads/Restraint — Pressure din bara de meniu superioară. 
Se introduce valoarea 5e7, în fereastra ce apare în partea stângă, ce reprezintă 50 - 10° 
N/m? = 50 MPa . Aplicarea condiţiilor de rezemare şi încărcare asupra piesei sunt redate 
în figura 13. 

Următoarea etapă o constituie discretizarea modelului piesei în elemente finite. 
Deoarece elementele finite sunt de tip solid, acestea au formă tetraedrică. Pentru 
realizarea discretizării se lansează comanda COSMOSWorks — Mesh — Apply 
Create... . Valorile ce apar în fereastra stângă se lasă nemodificate şi în acest moment se 
realizează discretizarea piesei în elemente finite, ilustrată în figura 14. 

După operațiunea de discretizare programul este gata de rularea efectivă a analizei cu 
elemente finite. Rularea se efectuează lansând comanda COSMOSWorks — Run din 
bara de meniu superioară. 
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In urma rulării analizei se trece la etapa vizualizării rezultatelor analizei, materializate 
prin afişarea grafică pe modelul piesei a stărilor de tensiune, respectiv de deformatie ce 
au luat naştere în urma solicitării piesei. 


Fig. 14. 
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Programul afişează câmpurile stărilor de tensiune orientate după axele sistemului de 
cordonate, precum şi tensiunile tangentiale, principale, sau tensiunile echivalente 
calculate conform teoriei a cincea de rezistență. Din punct de vedere al deformatiilor 
programul afişează câmpul de deplasări orientate după axele sistemului de cordonate, 
precum şi deplasările rezultante. În figura 15 este ilustrat câmpul tensiunilor echivalente 
calculate conform teoriei a cincea de rezistență iar în figura 16 câmpul deplasărilor 
rezultante. 
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